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THE EVOLUTIONARY DYNAMICS OF QUANTITATIVE CHARACTERS* 


P. NaRAIN 
Indian Agricultural Statistics Research Institute, New Delhi 110012 


INTRODUCTION 


The fact that evolution occurs is well recognised but what causes it to occur 
has been a matter of debate since the time Darwin® ascribed it to ‘Natural 
Selection’. With continued natural selection for several generations, individuals pos- 
sessing favourable characteristics have an advantage over those lacking them and as 
such they pass on such characteristics to their offspring for further modification and 
adaptation. Over the evolutionary time scale, this results in modified descendants of 
different ancestors living in geological times. However, for such a process to occur, it 
is necessary to have heritable variation and to understand factors governing it. Darwin 
unfortunately failed to discover the precise nature of such variation. It was only when 
Mendel’s laws were re-discovered in 1900, that the precise nature of such variation, in 
the form of gene mutation as raw material for evolution to act upon, was properly un- 
derstood. However, since the characteristics favoured during adaptive evolution are 
primarily quantitative in nature whereas mutation produce discontinuous variation, the 
compatibility between Darwinism and Mendelism was hotly debated for quite some time. 
Fisher’? paved the way for a reconciliation between the two by showing how the obse- 
rved correlation between relatives for such characteristics can be explained entirely in 
terms of the effects produced by Mendelian genes. Not only that, Fisher® further showed 
that the rate and direction of evolution are primarily determined by natural selection. 
According to his FundamentalT heorem of Natural Selection, the rate of increase in the 
average fitness of a population is proportional to additive genetic variance in fitness. 
Fisher*® also investigated, for the first time, the problem of the maintenance of genetic 
variability in natural populations for quantitative characters and gave a model in which 
fitness was assumed to decrease in proportion to the squared deviation from the opti- 
mum. In this he initiated investigation on a problem which was not only followed by 
several workers subsequently but is still being vigorously pursued since the mechanism 
for the maintenance of such variability in natural populations is far from being well 


understood. 


Subsequent to Fisher’s work, Wright®*”*!, Haldane’*, Robertson**, Latter** and 
Bulmer?*? dealt with models based on multiple loci but each with only two segregating 
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alleles. Modern understanding of molecular genetics, however, indicates that a gene is 
subdivisible into a very large number of variable nucleotide sites so that a model involv- 
ing segregation of infinitely many alleles at a given locus is more realistic than the early 
attempts. Crow and Kimura’ introduced, for the first time, an infinite allele model 
which was later studied by Kimura" in detail. In this study, mutation was considered 
as producing multiple alleles with varying phenotypic effects in a continuous manner at 
each of the several loci involved in the inheritance of quantitative characters. Later on, 
Latter?* adopted such a model in a discrete-time fashion. By extending this model to 
multiple loci, Lande'’~*1_ considered the effects of linkage and linkage disequilibrium. 
Fleming’® studied extensively the entire equilibrium structure of the multi-allele, multi- 
locus case, both for discrete as well as continuous time situations. 


The investigations of Latter?? and of Bulmer?’* based on diallelic loci, led to the 
conclusion that the equilibrium genetic variance was independent of the magnitude of 
the phenotypic effects produced by mutation but depended on the total mutation rate, 
summed over the loci, and the intensity of selection. On the other hand, in the infini- 
tely many allelle model of Kimura'‘, the equilibrium genetic variance depended on the 


underlying biological parameters and the distribution of allelic effects was found to be 
Gaussian approximately. 


The investigations of Turelli?’-** were based on an alternative approximation for 
the continuum-of-alleles model of Crow and Kimura® on the empirically motivated 
assumption that the effects of new mutations at a locus are usually much greater than 
the existing genetic variance at the locus. This led to the prediction of the same equi- 
librium genetic variance as those for diallelic loci, predicted by Latter?® and Bulmer?’. 


Subsequently, Nagylaki**, Gillespie!, Gimelfarb'? and Barton and Turelli! studied 
similar problems. 


In these studies (with the exception of Latter?’, Lande!’, Gillespie’, and Barton 
and Turelli'), the primary concern was to explain the magnitude of genetic variance 
within a population and little consideration was given to the genetic differentiation bet- 
ween populations or species. Moreover, none of these workers adopted a more realistic 
model of mutation involving discrete change of state. Chakraborty and Nei‘ however 
developed a new multation model called “discrete allelic effect’? to examine ie 
extent of genetic variation of a quantitative trait within a population as well as the 


same between two populations during the process of their genetic differentiation. But 
their study considered the forces of mutation and random d 


ignored. Narain and Chakraborty”® therefore, 


genetic saance within and between populations for quantitative characters determined 

by a few loci with major effects by using the discrete-time, discrete allelic-state model 

with mutation and selection in an infinitely large population. The selection was of th 

optimum type so that we could examine the change in variance under two situati : 

(a) a population that evolves from monomorphism (at an optimum phenotype) pe: 
, 


rift only and selection was 
studied the evolutionary changes of 
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(b) when an equilibrium population shifts to a new environment where the optimum 
phenotype is shifted by a few units of the phenotypic scale. The transient behaviour 
of the approach to equilibrium was studied in terms of changes in the means as well as 
the variances. The change in the interpopulational variance was examined when the 


equilibrium population splitted into two, where one of them moved to a new optimum 
environment. 


In this paper, we critically review some of the models referred to above for study- 
ing the dynamics of quantitative characters, both in terms of the means as well as_ the 
variances, under mutation-selection balance. In addition, we also study the genetic 
differentiation between populations or species in respect of quantitative characters. 


2. ALTERNATIVE APPROACHES 


The evolution of quantitative characters can be studied in either of the two ways. 
One way would be to define and study the underlying models at the level of phenotype 
avoiding any reference to gene frequencies. The Gaussian phenotypic analyses popula- 
rized by Lande'*’?° adopt this approach, on the assumption that genetic variances and 
covariances are known. Such an approach however, cannot predict the evolutionary 
dynamics of variance. 


The other way could be to define and study a genetic model assuming that a 
complete genetic analysis of the traits is possible. In the latter case, one has to start 
from the simplest situation of a single locus with two alleles and build over it the more 
complex systems of infinitely many alleles at a locus, several loci, linkage and 
epistatic effects etc. The results obtained from the simpler situations of one or two 
loci give an insight to the problem, particularly for characters such as skin pig- 
mentation in man which is believed to be controlled by a few loci say 5 to 6 with major 
effect. The rate of evolution, for such characters, is fairly rapid and modelling with few 
loci would be realistic. This latter approach would be mostly adopted in this paper. 


3. RELATIONSHIP BETWEEN PHENOTYPIC AND GENOTYPIC SELECTION 


The genetic properties of a given population are determined by gene and geno- 
typic frequencies and need to be connected to quantitative differences noticedin a 
metric character at the phenotypic level. We consider a population of individuals with 
k genotypes G,, G,,... G, with frequencies f,, fo, ..., fe and record the phenotypic mea- 
surements on each of the individual for the given character. Two individuals with the 
same genotype G; may then differ slightly in their measurements due to random effects 
ascribed to environmental differences. To account for it, we take, for the genotype Gi, 
a phenotypic distribution F (X;/G;) where X; is the random variable giving measure- 
ment on the individuals with genotype G:. The mean phenotype for this distribution 
is then 


g = J Fi f (Xi/Gi) aX. etl) 
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We thus have a sequence of k phenotypic distributions with means g), 22, ..., gx Corres- 
° . 2 
ponding to k genotypes but with a common environmental variance «o,. The average 


and variance of the phenotypic values of the character in the population are then, 


k 
M=2 fi gi 
‘-1 
k 2 
o, = eit (S58) th 9 po eee ---(2) 


the latter expression indicating that a, is the sum of the between-genotype variance 


(0% ) and the average within-genotypic variance (o7 ) due entirely to environmental 
effects. 


The change of the population mean resulting from the selection is brought about 
through the changes in the gene frequencies at the loci which influence the character 
under selection. But since the effects of the loci cannot be individually followed, the 
changes in the gene frequencies cannot, in practice, be ascertained unless we have some 
means of translating the phenotypic changes into genetic changes and vice-versa. The 
selection on the basis of the character X with a certain intensity induces selection among 
alleles at individual loci controlling the character. This is expressed in terms of a selec- 
tion function W (X) defined as relative selective value of an individual with measure- 
ment Y. Wecan then use the phenotypic distributions together with W (X) to determine 
the relative selective values W; conferred on the corresponding genotypes. This is simply 
the average fitness of individuals with a given genotype. These fitness values can be 
used, along with gene frequencies, to determine the gene frequency in the next gene- 
ration. This gene frequency information is then used to describe the mean and variance 
of the character in the next generation. The difference in the mean values between two 
successive generations gives the response to natural selection. 


We follow the method of Kimura and Crow’ to establish a general relation bet- 
ween the selection made at the overall Phenotypic level and the consequent selection 
induced at the genotypic level at individual loci. Let Yop be the optimum phenotypic 
value with maximum fitness. Consider a given locus with two alleles A and a 


, let the average phenotypic values of AA, Aa and aa indivi- 
duals be Tespectively X,,, X,. and Avot Measuring the character in units of phenotyic 
standard dlviation s,, we denote the density (before selection) and fitness functions b 

Ff (x) and w (x) respestively. If we take Xop as the origin, x = (X — Xop)/sp : 


Let m = (M — Xop)/o, and a1; = (Xi;—M) 
phenotypic value Y, j from the Population 
of the genotype with value X;,, then 


|p where a, is the deviation of the average 
meéan In gp units. Let w; y be the relative fitness 
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0) fe w (x) f (x — ary) dx. ...3) 


This is because the distribution of x in the sub-population of genotype with value Xi, 
is shifted by a;; and the density function in this sub-population before selection is 
f(x — a). This is a good approximation if the locus under consideration is contribut- 
ingonly a small portion of the total variance. Expanding f(x — a:,) in Taylor 
series about a;; = 0 and integrating, we get 


2 
ais 


Wiz = Bo— a1; By + > Be ...(4) 


where we neglect third order terms in a;,, treating it to be small, and 


hone i w (x) f (x) dx 


2 - 
B= Jf w(x)f' (x) dx et) 
—-oo 
"ie ° 
B= J w(x) f" (x) dx. 
— oO 
If the effect of substituting a for A is «, the average genotypic values of AA, Aa and 
aa can be expressed as a1, = — 2qa, a, = (1 — 2q)x and ao, = 2(1 —q)a respecti- 
vely. Then 
a; 


1 
Wir = Bo — a1, Bi + —> Be = Bo + 29 afr + 2q° a Be 


a 


1 = 2 2 : 
Wig = Bo—A2 Bi +5 B, = By — (1 — 29) 48, + ( 5 ) “ Be 
-..(6) 


2 
a 


2 
Wen = Bo — ny Bi + Br = Po — 2(1— ade Br + 21 — a) ot 
The mean fitness of the whole population is 
= (1 — gq)? wy + 2¢ (1 — 9) Wie + q? Woo 
= By +q(l — 927 Bs. AD) 
The gene frequency of a in the next generation is then 


gq =(q?W2t+ 9 il q) W1»)[® ---(8) 
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so that the change in gene frequency due to natural selection is 


Aq=q'-4 
= (g2 Wo + q (1 — g) Wiz — F@)/@ 
= q(1 —q)[— «Pi + a? Be (4 — Q)I/Bo = 59 (T ="q) an) 


neglecting terms involving «* and higher powers. This gives 


s = — « (Bi/B)) + 2? Be (4 — 9)]/Bo- ...(10) 


4. MUTATION-SELECTION BALANCE AT A DIALLELIC Locus 


Stabilising selection reduces genetic variability, but in natural populations this is 
countered in each generation by fresh variability generated by mutation. We have 
therefore, to determine how much genetic variability is likely to be maintained by the 
balance between these two forces. Following Bulmer’, we allow a mutation rate v from 
A; to A, as well as an equal mutation rate from A, to Ai. The expression for change. 
in gene frequency given by (9) then becomes 

Ag = q (1 — 4) [—« (B:/Bo) + 4 =? (Be/Bo) (1 — 2g)] + v (1 — 29). 
FES 


Let the fitness function be 
w (x) = exp[— } x*/o, ] ...(12) 


and suppose the mean is zero when g = 3 atall the loci controlling the trait. At 


equilibrium with g = q, the q's must be symmetrical about 1/2. 


Then 
(8,/Bo) = 0 
(8./B.) = — I/(6, +). : (13) 
This gives 
pa 
(1— 2g) [— 4a? q (1 — g)/(o6, + 52) +] = 0. ...(14) 


Thus q is eitherl/2 or it satisfies the quadratic equation 
A A ae 
a’q (1 — q) = 2v(o, + o% ). ...(15) 


Suppose now that there is an even number of loci (2k) and that (15) holds at all loci 
half of the gene frequencies being at the smaller and half at the larger root. Since 
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o *= 2ka? q(1 — q), we have 


> 


a” = 4ky(o%, + 3? (1 — 8kv) 
=4ky ‘Gs +o.) 


= 2kv/s. ...(16) 


As an illustration, suppose = 10c%. (weak selection) and that 4kv = 0.01 implying 


250 loci each with a mutation rate of 10-5, the expressed heritability under these rather 
favourable situations is only 0.10. The important point to note is that the equilibrium 
genetic variance is independent of the allelic effect. It depends only on the mutation rate 
and the intensity of selection. 


5. BALANCE IN AN INFINITELY MANY ALLELE MODEL 
The mutation model considered by Kimura’4 allowed for an infinitely many 
alleles at each locus. In this model, when a mutation occurs, it changes the value of 
the contribution of that allele by a small amount & from x to (x + &) where & is a 
random variable with density function f (&) which is symmetrical about zero with vari- 
ance o*. A continuous time model is further assumed so that the fitness of an indivi- 


dual with value x is — }x?/o” where the optimum value is set at zero. This is 
= p 


equivalent to the usual model of stabilising selection with fitness exp ( — } x?J0° ) in 


discrete time. 


Let p (x, t) be relative frequency of alleles with effect x at time ?. The rate of 
change of p (x, t) due to mutation is 


AA = — vp(x,th+y | p(x — & tf (8) dé 


6? p (x, t) ; 
ag op rer O (&) wltd) 


where v is mutation rate per generation, p (x — & f)is expanded around & = 0 by 
Taylor series and 


= yo? 


~ 1 
SE SE) ab = 0, | 


if Bf (E) da = 0, ...(18) 


ee ——+- — — —- 


is En f (E) dé = 0,n > 3. 


132 P. NARAIN 


The rate of change of p (x, t) due to selection is 


POD _ pyc, 1[— x2 + | x? p(x, 1) dx]o® 


or 


8—8 


= tp (x, t)[e; — x*]/o', ...(19) 
where 


ip xp (x, t)dx + 0 
—0o 


2 


x* p (x, t)dx > 


ss 


p (x, t) > p (x) 


as f — oo in the steady state. At equilibrium, therefore, we have, dp (x, t)/@t = 0 and 
2 mn ° 
$ vo? <2) — 4 p(x) (* — o° io, = 0. ...(20) 


This ordinary differential equation is satisfied when p(x) is the density function 
of a normal distribution with mean zero and variance given by 


3 => a) vor ...(21) 


Hence in the steady state, the balance between mutation of an infinitely many allele 
type and stabilising selection, in continuous time, produces a Gaussian distribution of 
allelic effects with zero mean and variance given by (21). The genetic value of the in- 





A 
dividual is therefore, distributed normally around zero mean with 20° ; 


6. BALANCE WITH A STEP-WiIsE DISCRETE MUTATION MopeL 
We consider a quantitative character controlled by 7 loci and assume that at each 
locus there is an infinite number of Possible allelic states. We assume that the pheno- 
typic effect of the alleles are discrete as shown in Fig. 1. 


In this figure, A; represents an allele occupying state i (any integer number from —co 
to oo) and having an allelic effect of ia. We assume that all allelic effects are additive 
with no dominance and no epistasis and that once A; mutates, it changes to allelic state 
Ai4, with probability. 

2m 


C= 4, = 
m—r 


ja forO << r<m ...(22) 


a, = 0, otherwise 
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-ma -2a -a 0 a 2a ma 
PHENOTYPIC EFFECT 


Fic. 1. Discrete allelic-state model used in this paper. In this model allele A; mutates to Ai,, 
with probability x, (= «_,), Allele Ai has a phenotypic effect of ai. 


where m is the possible number of mutational steps. The distribution is a shifted Bino- 
mial with mean zero and variance m/2. If v denotes the mutation rate, the absolute 
probability of such a mutation would be va,. Thus, an allele that mutates but has the 
same allelic effect as that of the original allele with probability va, so that in the con- 
ventional definition, the real mutation rate v’ would be given by v’ = (1 — ap) v. 
The per generation increment of the variance of allelic effect by mutation is then vma?/2. 


The selection operates on the total phenotypic value x and it is assumed to follow 


a normal distribution with mean p» and variance o7, . With no genotype-environment 


interaction, the distribution of x, among individuals of type A; A, with genotypic value 
a (i + j) follows a normal distribution with mean a(i +//) and the environmental 


variance c. with density function 


f (la i+ D) = oe exp[ — = ees | 23) 


The fitness function for the character value x is assumed to be Gaussian type as 


W (x) =Wmax. EXP |— 2 = Xen ...(24) 
20° 


where the character assumes the optimum fitness Wmax. at ¥ = Xopt. and Sy is the width 
of the function indicating the rate at which fitness declines with deviation of x from the 
optimum. Taking Wmax. = |, the mean fitness of the individuals with genotype A; A, 


would be 
wij = Jw (x) f (xJa (i + j)) dx 


= Oy V 2s exp [— s {a (i + j) —Xopt-}7] ni Zo) 


where s = 1/2 (o% + a ) indicates the strength of the selection at the group level. 


A large s, means weak selection of the stabilising type. 
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If x; (t) denotes the frequency of allele Ai in generation ¢ with allelic effect ai, an 
individual of genotype A; A, will have a mean reproductive fitness wi, = exp [—sa° 
(i + j)*] and thus, the change in gene frequency of Ai from generationt tot + 1 is 
given by 


dat)u(tt1=(l—-vt+ vao) & x1 (t) x; (t) exp [—sa® (i + j)*] 


+ y z Oy eo x, (t) [x14 (t) exp {—sa? (i + j + r)*} 


+ xi_, (t) exp {— sa* (i + j —r)*})) ...(26) 


where w, (t) is the mean fitness of individuals at the locus in the ‘th generation so 
adjusted as to make © x, (t) = 1. 
5] 


In general, this recurrence relationship does not yield any explicit solution. How- 
ever, for m = 1 it is possible to derive the equilibrium allele frequency profile by negle- 
cting powers of vy and s when the mean fitness ®@, (t) is approximated as 


Dya(th= l—s o° (t) ...(27) 
A 
where ay (t) is the total genotypic variance contributed by this locus at time? given 
by 
o,, (t) = a ~ Ex (t) x, (t)G +f). ..(28) 


The recurrence relation reduces to; 


(1) = > bn (t— + xa(t—- Di tu@—DIl—v—s 


OY = oe (eae (29) 
giving 


Axi (t) = ~» [xg — A OF xO 7 ne [ Sag? 


iceeacacegianmants — qrj2 Py 
2 


a Scalers wer ...(30) 
en ¢ population reaches equilibrium under the opposing pressures of mutation 
and selection, we have Ax; = 0. This gives 


A A 

x=(1 — SG) xo 

A A A 

*i4, — 2[1 —S(G— i*)x, + %_, =0,i > 1 


ix A 
Xi= Xi, 1 


...31) 


--—-~+--—-— 


non-zero integer. 
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The frequent allele frequencies are obtained as 
A A 

xX; = X (1 — SG) 

A 


A 
X, = xX [l — 2SG + 2S (1 —G) (1 — SG)] 
and 


xs = X9[l — 3SG + 6S (1 —G)(1 — SG) -..32) 
+ 2S {3 — SG (6 — G) + 2S(1 —G) 


x (4 — G) (1 — SG)}] 


W-———+,—* ————— 


A A 
ZUxx, (i +f)? = 0% [2a 
a) SA 


N]— 


where S = sa’/vyand G = 


In the general case of m-step mutational changes, the moments of the allelic effects 
as well as those of genotypic effects can be obtained analytically under optimum selec- 
tion. Denoting the kth moment of the distribution of allelic effects at a locus in the fth 


Co 
generation by M, (t) = 2 a*;*x; (t) and noting that for all i, x; (i) =x_: (t) at each 
+™-0O 


generation (since the optimum genotype is at origin), the recurrence relationship for the 
even order moments is given by 


m 2 
+ s[M, (t) Mex (t) — Mox+s (t)] 


Ma(teal) ={t =v + v( oi) (ie )e"] Mar () — M [Mox,2 (t) 


k 


+ (35 mod (2,) (4) oor 


1-0 
(43) 
so that the change of variance of allelic effects at generation rt, AM, (t) is given by 


mva? 


2 


At equilibrium, therefore, the fourth moment and the variance of allelic effects 
are related by 





AM, (t) = +s [M)(t) — Ma (t)]. ...(34) 


A 2 
Mp Mie oer : ...(35) 





Now, if pr (t) denotes the rth order moment of the genotypic effects at the locus 
(i. €., p(t) = BDal (i + J) x(t) x; (1), the variance and the fourth moment of the 
Le 


genotypic effect at a locus are related with those of allelic effects by 


Jy (t) = 2M, (t) and ps (t) = 2M, (t) + 6M; (1). 
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We thus obtain the change of variance of genotypic values at a locus in generation 
I, Ap, (t), as 


Apis (t) = mva* — 5 [us (0) — 2n2 ()). (36) 


A A 
If the equilibrium distribution of genotypic values is normal (i. e., uy = 3p, ), the ge- 


A A 
notypic variance ata locus inthe equilibrium population is given by a = pe 








= Vi mva*/s =a/ Ac’ |s where Ac., = mya?’ is the effect of mutational change on the 
genotypic value at a locus. This result is identical to that of Kimura** even though his 
model assumes a continuous distribution of allelic effects. 

Under the discrete allelic effect model, the departure from normality can be deter- 


A A A 
mined from the value of B, = y,/», at equilibrium. We thus obtain 
A ai ik 
B, = 2 + Ac™ /sa* 


A 
and hence the departure from normality as measured by y = 8, — 3 is given by 
A 
y = (Ac — so ss, aloe 


7. GENETIC DIFFERENTIATION BETWEEN POPULATIONS OR SPECIES 


Let us now consider the case where the populotion which is at equilibrium initi- 
ally with mean phenotype at origin and variance oF now shifts to a new environment 


where the optimum phenotyre is do» units away from the mean, obviously, under 
the effect of mutation and selection the equilibrium status of genotypic distribution 
will be immediately disturbed and gradually the distribution will shift towards the new 


optimum. To analyze the nature and rate of change of variability we must again con- 


sider the recurrence relationship of gene frequency changes. 


In this new environment, 
the fitness, 


Wij, Of an individual of genotype A; A; is given by 
Wy = exp [—s {ai + j) —d o,}). 


We then get, following similar deivations, 
4 +1) = wf (N—s [ps (t) + 2s (t) {ui (t) — do,}] ...(38) 


where yp’ (t) is the rth order moment (about origin) of the genotypic values in the ftth 


generation. Clearly, at equilibrium we thus have p’ 
1 


= ds,, the optimum genotypic 
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value and »; = 0. The equilibrium distribution being symmetric, the change in geno- 
typic variance is, similarly, given by 


Ao, (1) = mva* + s [20 (1) + ms (0) (dor) — 2, ()} — we 
(39) 


where 1; (¢), j44 (t) represent the third and fourth moments (about mean) of genotypic 
values at the ‘th generation. At equilibrium, since ,., is zero, we again have the same 
steady state genotypic variance given by 


Heres a ae aa i i ...(40) 


Thus, even if the optimum is shifted by a certain s. d. away from the original mean, as 
long as the intensity of selection (s) remains the same, the genotypic variance eventually 


returns to its original value although the genotypic distribution becomes now centered 
around the new optimum genotypic value. 





At the transitory stage however, it is difficult to assert analytically how the 
variance is altered. However, as we shall seein our numerical computations, at a 


transitory stage the genotypic variance first increases and eventually returns to its ori- 
ginal equilibrium value. 


It is apparent from the theoretical development given in the previous section that 
a quantitative study of genetic differentiation for metric traits between populations or 
species cannot be made analytically if optimal selection with stepwise mutation in an in- 
finite population is envisaged. A computer was therefore used to compute numerically 
the various quantities of interest by resorting to exact recurrence relations already dis- 
cussed in the previous sections, The mean and variability of the character both within 
as well as between populations were studied in the transient stage and at equilibrium. 
To start with, the initial population was considered as monomorphic at optimum and 
the behaviour of within population variance studied over time. After reaching the 
equilibrium, the optimum was shifted to a few units on the right and the transient be- 
haviour of the mean and variability (within as well as between) was studied. 


Since the mean is at the optimum and the optimum is set at zero, the population 
mean remains at zero unless there is a shift of the optimum genotype. However, the 
variance increases slowly from zero and attains, at equilibrium, a value determined 
solely by (v/s) as already obtained algebraically. This behaviour of within population 
variance as a function of time for different values of s between 0.004 to 0.040 for 

= ] and v = 0.0005 is presented in Table I. 


For intense selection, the variances during the transient stage as well as at equilibrium 
are lower as otherwise expected. Mutation creates variability while selection eliminates 
it so that for intense selection, its role is dominant. Also, the approach to equilibrium 
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TABLE I 
Within population variance at equilibrium when initial population is genetically homogeneous 
with only one type of individuals which is optimum phenotype for some selected values of 
selection coefficients (s). The mutation rate, v, is taken as 0.0005 andm = 1 


LL ——EEEEEe 
Ae 


Ss S94 
Tapia e in ij: 000A brary oy =” bok ian leit en 
0.008 6.03 X10~* 
0.010 4.85x1072 
0.020 2.44107? 
0.040 1.22X10~ 





is quicker for more intense selection as revealed by results of computer simulations 
which are not presented here. 


When we shift the optimum to four or six standard deviations away from 
the mean on the right and study the transient behaviour of the process as it approaches 
the same equilibrium, we notice some interesting results. In Table II, we present these 
results in terms of mean, variance, skewness and kurtosis of genotypic values in 
different generations for a quantitative character under centripetal selection when the 
new optimum is at six standard deviations away from the original optimum. We take 
v= 0.001, s = 2v, and m = 5. 


TABLE II 
Mean (X;,), variance (V.), skewness (y\z) and kurtosis (y2t) of genotypic values in different gene- 
rations for a quantitative trait under centripetal selection (v =0.001, s= 2v and m = 5) with new 
optimum 6 standard deviations away from original optimum 


eS 


t 


0 50 100 500 1000 re) 
pg ee ee 
Mean( 11) 0 0.378 0.930 3.443 3.811 4.000 
Variance (V;) 0.456 0.728 1.296 0.815 0.618 0.456 
Skewnessy; (¥11) 0 0.997 0.684 0.119 —0.001 0 
Kurtosis (‘y2:) 2.741 1.886 0.433 0.764 1.563 2.741 


eee 


N. B.: Kurtosis (y:+) is measured by y2: = (Hae/V, y-3, 


The mean increases from zero to four at equilibrium but the increase is more Tapid 
for a higher value of m. The variance on the other hand increases, attains a transitory 
maximum and decreases back to the original steady state value which is related to the 
fourth moment of the genotypic values as given by (40). When s is very large compared 
to v, more than one transitory maxima are produced. In Fig. 2, we present such a 
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TIME (t) 
Fic. 2. Intra-population variance (Vi) under the joint effects of mutation and centripetal 
selection as a function of time /, with m = 1, v = 0.001, s = 0.002, 0.004 and 0.008 


when the new optimum is at 6 standard deviations away from the original optimum at 
the origin. Initial population is at steady-state. 


behaviour of intra-population variance for s = 0.002, 0.004 and 0.008 when v = 0.001, 
m = 1 and when the new optimum is at six standard deviations away from the original 
optimum at the origin. 


For a higher value of m, the variance attains a considerably higher peak as well 
as somewhat earlier than when m = 1. The most interesting feature is regarding the 
skewness of the distribution of genotypic values. Initially, this distribution is symmet- 
rical but as we advance in time, its symmetry is disturbed. It gets skewed initially and 
then slowly the skewness decreases, changes sign and finally the distribution becomes 
again symmetrical at equilibrium. The kurtosis of the distribution also behaves ina 
similar fashion. Starting from a value very near to three initially, it declines to a value 
less than half but increases thereafter and restores the initial value at the equilibrium. 
Compared to m = 1, the distribution for m = 5 becomes more leptokurtic. A typical 
allelic distribution corresponding tov = 0.001, 5s = 2vandm = 5, depicting these 
features, is shown in Fig. 3. 


Of special importance in these studies is the genetic differentiation between popu- 
lations built up over a period of time when in one population the same optimum holds 
but in the other it has shifted a certain distance away from the mean. Chakraborty and 
Nei‘ used the ratio of between population (B,) to within population (V;) variability as 
an index for determining the evolutionary forces under which the character changes 
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Fic. 3. Allele frequency distribution under joint effects of mutation and centripetal selection at 
different generations (t = 0, 60, 150, 1000 and ec) with y = 0.001, s = 2vand m= 5 
when the new optimum is 6 standard deviations away from the original optimum at the 
origin. Initial population is at steady-state. 


over time. Their studies involving drift only showed that this ratio (B,/V;) increases 
linearly with time. With centripetal selection in an infinitely large population . we find 
that this behaviour changes considerably and it is no longer a monotone function of 
time. 


We have already discussed the behaviour of V, which increases slowly from initial 
equilibrium value, attains a maximum and then decreases back to the same value at equi- 
librium, But when we consider between population variability, B,, it is found that it in- 
creases slowly initially and then almost linearly until it approaches a plateau at equilibri- 
um. The ratio (B,/V,) almost mimics the behaviour of B, at least in the initia] stages but 
it attains a much higher value. This is obvious because V; decreases while B; increases 
as equilibrium is reached. In the initial transient stage, however, B, and the ratio are 
almost the same because V, has been increasing and reaching a maximum. After this 
stage, at which maximum JV, occurs, the quantities B, and (B,/V,) diverge apart, increasing 
with time by different magnitudes. To illustrate the qualitative nature of the changes 
in V;, B, and (B&/V,), under joint effects of mutation and centripetal selection, 
Fig. 4 presents the numerical results for vy = 0,001, s = 2y and m = 5 when the new 


optimum is taken to be approximately six standard deviations away from the optimum 
in the other population. 


It is seen from this figure thatjthe rate of approach to equilibrium variance within 
population (0.456) is faster as compared to that of the mean Senotypic value (4.0) since 
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VARIANCE AND VARIANCE RATIO 





TIME (t) 

Fic. 4. Intra- (V;) and inter-populational (B;) variance components under the joint effects of 
mutation and centripetal selection as a function of time of divergence of two populations 
one of which has optimum genotye six standard deviations away from the other (which 
remains at a steady state by mutation-selection balance). Time (t) is measured in 
units of generation. Initial population is at steady state. The parameter values are 
vy = 0.001, s = 2vandm = 5S. 


the between-population variance (B;) attains its equilibrium value (4.0) at a later 
time than V;. The variance ratio (B,/V,) reaches a steady state value (8.80) since the 
process of genetic differentiation stops once the diverging population reaches its steady- 
state genotypic distribution around its new optimum value. We thus see that (B,/V;) 
asymptotes and does not therefore increase linearly with time as in Chakraborty and 
Nei‘ under mutation-drift balance. The behaviour of (2,/V;) as a function of time of 
divergence can therefore be taken as a criterion for determining whether selective forces 


are operating or not. 
8. DISCUSSION 
Kimura™ showed, for infinitely large populations, and assuming a continuous time 


process that the distribution of allelic effects tends to be normal at equilibrium bet- 
ween selection and mutational forces and that the mean and variance of the equilibrium 
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distribution are determined by the amounts of increase in mean and variance of the 
genotypic value per gene per generation as well as by the intensity of fitness function. 
The discrete allelic-state model with the assumption of a discrete-time process, consi- 
dered in this paper, has revealed behaviour similar to those of Kimura‘ as it should, 
since binomial distribution of allelic effects should tend to normal distribution as we go 
from discrete to continuous case. On the other hand, the diallelic model of Latter??, re- 
analysed by Bulmer®’® who used different approximations but arrived at the same con- 
clusions as those of Latter**, indicated different results. Their predictions for equilibrium 
genetic variance differ qualitatively from those of Kimura-Lande-Fleming-Narain & 
Chakraborty. 


Based on “‘House-of-cards” approximation of Kingman?®, Turelli2? presented a 
new asymptotic analysis of Kimura’s model to show that the qualitatively different pre- 
dictions about equilibrium genetic variance are not due to the number of alleles assumed 
per locus. Instead, such different results are attributable to assumptions concerning the 
relative magnitudes of per locus mutation rates, the phenotypic effects of mutation 
and intensity of selection. He then analysed a model with tri-allelic loci which 
allows among-locus variation in mutation rates and allelic effects. At the single 
locus level, such a model is a particular case of the discrete-time, discrete allelic-state 
mutation model analysed in the present paper and leads to the same conclusions as in 
Turelli*’ in so far as the equilibrium genetic variance is concerned. 


As already mentioned in the Introduction, very few studies on the problem of 
genetic differentiation between population or species have appeared, particularly for 
the situation when one of the daughter populations has a shifted optimum. Such cases 
have biological relevance as for instance in skin pigmentation for a small group of 
Caucasian race with fair skin who moved out of Central Asia around 3000 years ago 
and settled in southern parts of America with plenty of sunlight. The optimum pheno- 
type for skin pigmentation must therefore, have shifted by several standard deviations 
away from the original optimum. This introduces differentiation between populations 
and it is of interest to study the transient properties of such a process. This has been 
done in this paper by Studying the distribution of allelic frequency asa function of 
time. Significant changes in statistical Properties of the distribution such as mean, 
variance, skewness and kurtosis have been noticed. In particular, algebraic expressions 
for changes in mean and variance reveal interesting results. When we consider optimum 
at ; standard deviations units away from the origin, the mean change in the genotypic 
value of the character at the f-th generation, as given by eqn. (38), depends on s, d, Sp, 


Ms (t — 1), wo (t — 1) and Hy (t — 1). It does not depend on the mutation rate, v. But 
change in the genotypic variance per generation, as given by eqn. (39), clearly indi- 
Cates that it is affected by the mutational component. Initially, the population is in 


equilibrium and symmetrical. If we further assume that it has normal kurtosis, we have 


#1 (0) = 0, us (0) = 0, py (0) = 3p, (0). If we also disregard the increase in varia- 
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tion due to mutation, the expression (39) reduces to the result given in Latter?*. The 
genotypic variance is reduced and the reduction depends on the intensity of selection 
and the heritability. 


The transient behaviour of the mean noticed in this investigation can be useful in 
giving some idea about the divergence time; at least the time by which the mean gets 
to half of the total change in the mean, i. €., 1/2 (doy), which is known. The computer 
results show that intense selection speeds up the divergence time but by increasing the 
number of mutational steps, this time is shortened, though not very appreciably. Thus, 
for v = 0.001, s = 0.008 and m = 1, this time is of the order of 120 generations, With 
s = 0,002 and m = | and 5, the times are respectively 300 and 170 generations. 


The main focus of this study is on the genetic differentation between populations 
which gets built up over time when the population, after reaching equilibrium by mut- 
ation-selection balance, splits up into two in one of which the same optimum holds but 
in the other it shifts a few standard deviations away from the mean. The ratio between 
versus within population variance (B,/V;) as a function of time then provides with a 
possible mean of ascertaining the role of adaptive changes under which the character 
changes over time. With selection, this ratio necessarily changes non-linearly with time. 


SUMMARY 


Studies on the maintenance of genetic variability for a quantitative trait due to a 
balance between stabilising selection and mutation in natural populations have been re- 
viewed. Several models, both in terms of the dynamics of the means as well as the 
variances have been discussed. In particular, using a new and more general genetic model 
called the discrete-allelic state model and assuming discrete-time process, the evolutio- 
nary changes of genetic variation of quantitative characters, controlled by a few loci, 
within and between populations during the process of genetic differentiation of popula- 
tions or species, are studied under the effects of mutation and centripetal selection in 
infinitely large populations. While in a finite population and ignoring selection, the 
rate of change of additive genetic variance depends on mutation and effective population 
size, traits under optimal selection in infinitely large populations go through the dyna- 
mics of a rather complicated form depending on the relative intensities of selection aud 
mutation. When a population, which has reached steady-state by mutation-selection 
balance, splits into two, in one of which the same optimum genotype holds but in the 
other the optimum shifts a few standard deviations away from the original optimum, 
the corresponding daughter population starts differentiating from its sister population 
by favouring certain class of mutant alleles and discarding others which were originally 
favoured. During this process of turn over of genes, both the intra- and inter-popula- 
tional variances undergo a complicated change, and the ratio of the latter to the former 
is a non-linear function of time of divergence. This pattern is qualitatively very different 
from the case when selection is absent. The intra-population distribution of genotypic 
values, during this transition, is shown to deviate considerably from normality. 
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Let /(n, k) denote the number of ways of selecting k balls from a balls 
arranged in a line (called an 7-line) with no two adjacent balls, from the k 
selected balls, being unit separation. It is shown that /(n, k) satisfies the re- 
currence relation 


lan,k)= Il(n-1,) +1(2-1,k -1)—Il(n-2,k —-1I)+1l(n—3,k -1) 


and an explicit form for /(n, k) is obtained. 


1. INTRODUCTION 


Kaplansky' proved that, 


ee aes) o<ke i! 
f (nk) = k 


0 otherwise, 





where f (n, 4) denotes the number of ways of selecting k balls from n balls arranged in 
a line without any two selected balls being consecutive, (see also Riordan’®, Ryser®). 


Recently Konvalina® proved that 


f [k/2] A k —2i 
—k-2i\... — 
g(n,k)= | (" k — 2i )ifms2@-0 
1 g=0 
ro ifn 2 (ke 1) 


where [k/2] is the greatest integer less than or equal to k/2 and g(n, k) denotes the 
number of ways of selecting k balls from balls arranged in a line (called an n-line) 
without any two selected balls being uni-separate (i. e. being separated by exactly one 
ball which can be either selected or not). 


A related problem is to determine the number of ways of selecting & balls from 
an n-line with no two adjacent balls, from the k selected balls, being uni-separate. Let 
1 (n, k) denote the number of ways of selecting k balls from an n-line with no two adja- 
cent balls, from the k selected balls, being uni-separate. 
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In this paper we derive a recurrence relation for /(n, kK) and then give a method 
similar to that used by Riordan‘ to determine / (a, k). 


2. THe MAIN RESULT 
Lemma 2.1—Let / (n, 0) = 1, then / (n ,k) satisfies the recurrence relation 


1 (n,k) = 1(n — 1,k) + 1(n —1, kK—-1) — 1(n —2,k-1) + 1 (n—-3, &-1) 
(2:1) 


with the boundary conditions / (n, 1) = nand/(1, k) = 0Oifk> 1. 


Proor: To prove the recurrence relation (2.1), let / (n, k) be defined as before. 
Then, the corresponding selections either contain the first ball or they do not. If they 
do not, then they are enumerated by / (n — 1, k). If they do, then the selections either 
contain the second ball or they do not. If they do not, then they cannot contain the 
third ball and, hence, are enumerated by / (n — 3, k — 1). If the selections contain the 
first and second balls, then the selections either contain the third ball or they do not. 
If they do not they cannot contain the fourth and, hence are enumerated by 
I(n — 4,k — 2). If the selections contain the first, the second and the third balls, 
then the selections either contain the fourth or they do not, and so on. Hence, 
I (n, k) satisfies the recurrence 


I(n,k) =1(n—1,k) +1 (n — 9, Kise 1) “0 (8 2 4) ee Oo) 


= I(n — lk)+](n—3,k —- 1) + L (n, k) 
where 


k 
L(nk)= S1(n—2—i,k — 3). 


Then 
I(nm-—1,k— N=I(n-2k — 1) +1 (n -—4,k — 2) 
+l -—-5,k—3) +... 
=I(n-—2,k — 1) + L (a, k) 
thus the recurrence relation (2.1) is obtained. 
Theorem 2.1— 


a 


(nk) = > (iste) (apenas ~'YVitk<n 


i=0 


where 





A= min( &—1,[ 2 )) 
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and / (n, k) = 0 otherwise. 


PRooF : First we find / (n, k) for k = n,n — | andn — 2. Itis clear that 1 (n,n) 
= 1, from the boundary conditions and using the recurrence (2.1), we obtain 


I(a,n—1)=12— 1n—1)+1(n— l,m — 2) — I(n — 2, — 2) 
+ I(n — 3,n — 2) 
=I(n— 1,n—2) =... = 1(2,1) = 2 
I(n,n—2)=I1(n—1,n—2)+1]Qn— l,n — 3) —I(n —2,n — 3) 
+ 1(n — 3,n-- 3) 
=I(n -—1,n—3)+1 
= 13,1) +n—3 
=3+n—3=—n,. 
It remains to find the explicit form for / (n, k) if k <n — 3. From the boundary con- 
ditions / (n, 0) = 1, /(n, 1) =n, and using the recurrence (2.1) we have 
I(n, 2) =1(n — 1,2) + 17 (n — 1,1) —1(# — 2,1) 4+1(n — 3, 1) 
= I(n — 1,2) +(n — 2) 
= 1(2,2)+(n—2)+(n - 3)+...4+1 


(n — 1) (n — 2) __ (n—2) (n—3) 
ae el + 1 eels cer aeett 


mr (Vota) estes) 


(n—1) 


Similarly 


Peay (ey 2) (9 — 22) h(n "3; 2) 
mel eT Set ae (ieee) Ft famed + (n —3) 


La+m—4 + ("> °)+e-3) 


= ](n 


=I1(n 


13) 4 fe F *) (n— 2) 


which entails 


ren 3y—= ("7 7) +2("2 °) + ("5") 
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By the mathematical induction, we get 


d 
k—1\(n—k+1-i). Ee Pe 
inb= > ( j )( i+] ) it <nA = mink l, 


(¢ )): 


It is clear that / (n, k) = O if k > n, this completes the proof of the theorem. 





Finally a table is given for the numbers / (n, k) where n = 0 (1) 12. 
The numbers / (n, k) 





kIn 0 | 2 3 4 = 6 cj 8 ee) |e? Res F 
0 ] | ] ] 1 ] ] ] 1 1 ] 1 1 
1 1 2 3 4 5 6 7 8 9110 a i- St? 
2 l 2 4 fe A 162 22) 29 a aos 
3 ] 2 5. 10. IS | 30." 47 = FOS ise ae 
4 | 2 Gy 13. 326%) AS 839146 io 

5 | 2 7 169 3355 eee 
6 l 2 8 19 45 96 192 
7 I 2 22 56 126 
8 ] 4, 0107225 68 
9 I 2° 1 aes 
10 1 2 eis 
k, 2 
12 1 

a ee 
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It is shown that a real linear metric with Pseudo strict convexity and metric 
convexity induces a strictly convex norm. The case of a complex linear metric 
is also completly analysed. 


Ahuja et al.' and Sastry and Naidu*' studied extensively the mutual implications 
among various convexity conditions in metric linear spaces. In that connection, a 
natural question that arises is whether a metric linear space with pseudo strict conve- 
xity and metric convexity turns out to bea normed linear space. We answer this 
question in the affirmative in the real case and analyse the complex case. The remark 
on page 554 of Assad and Kirk’ plays a key role in the proofs 


Let (X, d) be a metric space. We say that d is metrically convex (m . c) if given 
x # y, there exists z(~ x and y) such that d(x, z) + d (z, y) = d(x, y). 


Suppose (X, d) is a metric space with metric convexity. Fix x, yin X. Write 
S = {z € X/d (x, z) + d (z, y) = d (x, y)}. 


Clearly Sis a closed set. Foru, vin S, define u < v if d (x, vy) = d(x, u) + d (u, v). 
Then < isa partial order on S._ Let A be a maximal totally ordered subset of S. 


Then A has the following properties : 
(I) x,» € A and A is closed. 
(If) Ifu, v € A andu < », then there exists w © A such thatu < w < ». 
(III) Suppose (X, d) is complete. Define 9: A > R as 9 (z) = d (x, z). 


Then gis a strictly increasing continuous map of A onto [0, d(x, y)]. 
Further 9 is strictly increasing and continuous. Thus 4 is a perfect set. 
Infact, ifz € A and x < zy then there exists a strictly increasing 
sequence in A converging to z; if x < z < y, then there exists a strictly 
decreasing sequence in A converging to z. 


150 K. P. R. SASTRY ET AL. 


As a consequence of this, we have 


(IV) (Assad and Kirk’). Suppose (X, d) is complete with metric convexity. If K 
is a closed subset of X,x € K and y & K, then there exists z € OK 
(boundary of K) such that d (x, z) + d (2; y) =. d.(x,y). 


The following example shows that the completion of a metrically convex metric 
space need not be metrically convex. It also shows that (IIT) and ([V) above may fail 
if the metric is not complete. 


Example—Let X = [0, 1/6] U (1/3, 4] U (2/3, 1] with the usual metric on the 
real line. 


Definition—Let (X, d) be a metric linear space. We say that 


(i) dis pseudo strict convex (p.s.c) if x #0, y 4 0, d(x + y, 0) = d(x, 0) 
+ d(y, 0) imply that y = tx for some t > 0. 


(ii) d is ball convex (b.c) if d(x, 0) = r = d (y, 0) implies that d (4 (x + y), 0) 
< r (that is, balls are convex). 


(iii) d is strictly convex (s.c) if r> 0, d(x,0)<r,d(y,0) <r imply that 
d (4 (x + y), 0) < r (that is, the balls are convex and do not contain line 
segments on the surface). 


Lemma 1—Suppose d is a linear metric on R, the real line. Then d defines a norm 
if it is metrically convex. 


Proor : Suppose (R,d) is ‘metrically convex. Suppose 0<t<s._ Let 
K =[-—s,t]. Then K is closed in(R,d),0 € K ands & K. Hence, by (IV), there 
exists z € 0K such that d(0, z) + d(z,s)=d(0,s). But —s and ¢ are the only 
points of dK. Hence z= tf, so that d(0,1)+ d(t,s) =d(0,s). Thus, d(0, x) 
+ d(0,y) = d (0, x + y) for all non-negative x,y. From this we conclude that d (0,t) 
= td(0, 1) for allt > 0. Since d(0, t) = d(0, —12), this leads to the conclusion. 


Lemma 2—Suppese (X, d) isa real metric linear space. Then d defines a norm 
on X if and only if every one dimensional subspace of XY has metric convexity. 


PRooF: Applying Lemma ] to one dimensional subspaces, the result follows. 


Lemma 3—Let (X, d) be a real or complex metric linear space. Suppose that 


(X, d) has metric convexity and pseudo strict convexity. Then every one dimensional 
subspace of (X, d) has metric convexity. 


Proor : Let Y be a one dimensional subspace of (X,d). Suppose x, »y € Y and 
x fy. Then, by metric convexity of (X,d), there exists z (~ x and y) such that 
d (x, z) + d(z, y) = d(x, vy). Hence d(x — 2,0) + d(z — y,0) = d(x — y, 0) 
Since (X, d) has (p.s.c), it follows that z — y = 1 (x — z) for somer > 0 Lanesnoatiis 
zisinY. Thus ¥Y has metric covexity. 
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Theorem 4—Let (X, d) be a real metric linear space. Then d defines a strictly 
convex norm on X if and only if (X,d) has metric convexity and pseudo strict 
convexity. 


Proor : If (X, d) has (m.c) and (p.s.c.), then, by Lemma 2 and Lemma 3, d 
defines a norm. Since, ina normed linear space, (s.c) and (p.s.c) are equivalent, it 
follows that d defines a strictly convex norm on X._ The converse is evident. 


The case of a complex linear metric is different. The following example shows 
that in a one dimensional complex metric linear space with (m.c) and (p.s.c.), the 
metric need not define anorm. In what follows, C is the field of complex numbers, 
regarded as one dimensional complex linear space. 


Example 2—Define the linear metric d on the one dimensional complex linear 
space C as 


d ((0, 0), (a, b)) = (Ja/® + /b/’)/3 for alla, bin R. 
Then (C, d) has (m, c), (p.s.c) and (s.c) but d does not define a norm on C. 


The following lemma characterizes norm-defining linear metrics on one dimen- 
sional complex linear spaces. 


Lemma 5—Suppose d is a linear metric on the one dimensional complex linear 
space C. Thenddefines a normonC if and only ifd is metrically convex and 
d (z, 0) = d (/z/, 0) for all z in C, where /z/ is the usul absolute value of z (that is, 
metric is rotation invariant). 


Proor : Suppose d is metrically convex and rotationinvariant. LetO<t<s. 
Write K = {z © C//z] < t}. Then Kis closedin (C,d),0 € Kands & K, so that 
by (IV), there exists w € OK such that d (0, w) + d(w, s) = d(0,5). Since w € @K, 
we have /w/ = t. Since d (0, w) = d (0,/w/) = d (0, t), we have 


d (0, t) + d(w,s) = d (0, s). (1) 

Thus d (0, t) < d (0, s) whenever 0 <t<-s. 

Since s — t = /s/ — /w/ & /s — wi, we have 

d(0,s —t)<d(0,/s — w/) = d (0, s— w)so that d(t, s) < d (uv, s). 

This, together with (1), gives d (0, t) + d(t, s) = d (0, s). 

Consequently, d (0, a) + 4(0, 5) = d (0, a + b) for all non-negative a, 5. 

Now, using the continuity of d, we conclude that 

d(0, a) = ad (0, 1) for all non-negative a. 

Finally, for a, z € C, ; 
d (0, az) = d(0,Jaz/) = Ja]|z[d (0, 1) = Jal d (, |z/) = |al d (0, z) 
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showing that d defines a norm on C. 
The converse is evident. 


The following theorem characterizes norm-defining linear metrics on complex 
linear spaces. 


Theorem 6—Let (X, d) be a complex metric linear space. Then d defines a norm 
on X if and only if (i) every one dimensional subspace of (X, d) has metric convexity 
and (ii) d (ax, 0) = d (/a/x, 0) for alla € C and x € X. (condition (ii) is superfluous 
in real metric linear spaces). 


Proor : Suppose (i) and (ii) hold. Then Lemma 5 shows that d defines a norm 
on every one dimensional subspace of X. That is, d defines a normon X. The 
converse is evident. 


The following theorem is a consequence of Lemma 3 and Theorem 6. 


Theorem 7—Let (X,d) be a complex metric linear space. Then d defines a 
strictly convex norm on X if and only if (i) (X, d) has (m.c) and (p.s.c) and 


(ii) d (ax, 0) = d (/a/x, 0) for alla € Candx € X. 


Lemma 8—Let (X, d) be a complete metric linear space with metric convexity 
and ball convexity. Then d (x, z) + d(z, y) = d(x, y) for all x, y in XY, where z =} 
(x + y). In particular, every one dimensional subspace of (X, d) has metric convexity. 


Proor: Let x,y € X. Write r = d(x, y). By (III), there exists w € A 
such that d (x, w) = r, so that d (w,y) =r. Now, 


d(x,x+y —w)=d(w,y)=r=d(x,w) =d(y,x+y— w). 


Thus w and x + y — w are in both the balls B(x, r), and B(y, r) where B (x, r) 


= {ufd (x, u)€r} and B(y, r) has similar meaning. Since (X,d) has ball convexity, 
it follows that 


z=4(w+ 2z—we B(x,r)M B(y, r). 


Hence, d (x, y) < d (x,z)+d(z,y\or+r=d (x, »), which gives the result. 


As a consequence of the above lemma, we h 


ahd ave the following characterization of 
norm-defining linear metrics. 


Theorem 9—Let (X, d) be a com 


plete metric linear space. Then d defines a 
norm on X if and only if 


(i) (X, d) has metric convexity and ball convexity and 
(ii) d(ax, 0) =d (/a/x, 0) whenever a is a scalar and x € Y. 


(condition (ii) is superfluous in the real case). 
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The proof follows from Lemma 2, Theorem 6 and Lemma 8. 


The following example shows that a metric linear space need not have ball 
convexity, even if every one dimensional subspace is strictly convex. 


Example 3—For (x, y) € R?, define 


_ Jmax {/x/, /yf} if /x) + ly] <1 
4 (x; y), (, OY) ee fIx/, [oP (+ [xl -+IvDI2/x] +2)y)) 


otherwise. Then (R’, d) is areal metric linear space which does not have any of (mc), 


(b.c) and (p.s.c). But every one dimesional subspace is unbounded and strictly 
convex. 


The following is an example of a metric linear space without ball convexity, in 
which every one dimensional subspace has ball convexity. 


Example 4—For (x, v) € R?’, define 


[x/+]y] if /x/+/y] <4 
d ((x, y), (0, 0)) = 23 if /x/+/y/> 4 and /y/<3 or if O</y/—43<[x/ 
max {/y/—/x/, 1} if [y[—|x/> 3. 


Then (R?, d) has neither ball convexity nor (p.s.c), but every one dimensional subspace 
has ball convexity. While some of the one dimensional subspaces are bounded, some 
others are unbounded. 


Problems : (1) Does a metric linear space with (m.c) have (b.c) ? 


(2) Is the completion of a metric linear space with (m.c) and (b.c), metrically 
convex ? 
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Let X be acompact Hausdorff space. Let +:X-—-X bea homeomorhpism 
such that t? is the identity on X. Let C (X) be the Banach algebra of all 
complex-Valued continuous functions on X¥. Let C(X,7) ={fEC(X): 
f(:(x)) = f (x) for all x € X}. In this paper the following characterization 
of C (X, ) is proved : 

If A is a real function algebra on (X, =) and if Re A is uniformly closed, then 
(Alex CCX"): 


1. INTRODUCTION 


Let X be a compact Hausdorff space, C(X) the Banach algebra of continuous 
complex-valued functions on X with the supremum norm and A, a function algebra on 
X (that is a uniformly closed subalgebra of C (X) containing constants and separating 
points of X). The classical Stone-Weierstrass theorem gives a condition (namely that 


f © Aimplies f € A) which is necessary and sufficient for Ato be the whole of 
(sie af 


Since the appearance of the Stone-Weierstrass theorem, several other conditions 
have been discussed. A systematic account of these can be found in Burckel!. Many 
of these conditions are in terms of Re A, the space of real parts of functions in A. 
One of these characterizations is given by the Hoffman-Wermer theorem which states 
that if Re A is uniformly closed, then A = C (X). (Hoffman and Wermer’). 


In the present note, we prove an analogue of the Hoffman-Wermer theorem for 
a real function algebra. This states that if Yis a compact Hausdorff space, + is an 


involutoric homeomorphism on_X, A is a real function algebra on (X,7) (see Definition 
2.1), and Re A is uniformly closed, then 


A= C(X,+) = {f € C (X): f(s (x) = F(x) forallx E xv 
(Theorem 3.1), In proving this theorem, we have made u 
theorem proved earlier by us* 

subsets. The next section con 
sets and antisymmetric sets. 


se of an analogue of Bishop’s 
which deals with the partition of X into antisymmetric 
tains the basic definitions and some properties of peak 
These properties are used in the last section to prove the 


ANALOGUE OF HOFFMAN-WERMER THEOREM 155 


main theorem. This main theorem is then used to prove the following: Let A bea 
uniformly closed real subalgebra of C (X), containing real constants and separating 
points of X. If Re A is closed, then there exists a closed subset Z of X such that 


A={fEC(X): f| Z is real} (Theorem 3.7). 


This can be regarded as a stronger version of the classical Hoffman-Wermer Theorem 
(Remark 3.8). 


We have undertaken a systematic study of the conditions on a real function 
algebra A on (X, t) which force A to be the whole of C(X, +). The present note 


contains one result in this direction. For other results along similar lines, see 
Kulkarni and Srinivasan’. 


2. ANTI SYMMETRIC SETS AND PEAK SBETs 


Let X be a compact-Hausdorff space and C(X), the Banach algebra of all 
complex valued continuous functions on X with the supremum norm. 


Let t: ¥ ~ X be a homeomorphism such that t? = + © t is the identity on X. 
Let 


C(x, =) ={7e€ C (X) -Sf (= (x) = f () for all x € X}. 


Then C (X, +) is a real communtative Banach algebra with the identity 1. Also it is 
easy to see that C (X, t) separates the points of X, that is for any x1 # x, in X there 
is f in C (X, +) such that f (x,) # f (%)). 

Definition 2.1°—Let A be a subalgebra of C (X, 1). Then A is said to bea real 
function algebra on (X, +) if (i) A is uniformly closed; (ii) A separates points of X; 
(iii) A contains (real) constants. 


Note that, if f € A and / is constant, then f must be real. 


Definition 2.2°—Let A be a real function algebra on (X, 7). A nonempty subset 
K of Xis called a set of anti-symmetry on A or A-antisymmetric if (i) f € A and 
f | K is real implies f | K is constant; and (ii) f € A and f | Kis purely imaginary im- 
plies f | K is constant. 

For the properties of A-antisymmetric sets we refer to the earlier paper of the 
authors cited above. 


Definition 2.3—Let A be a real function algebra on (X,+). Then 4 is said to 
be an anti-symmetric algebra if X is an A-antisymmetric set. 


Definition 2.4—Let A be a real function algebra on (X, +). Let S be a non-empty 
subset of X then S is said to be A-peak set if there exists f € A such that 


S = {x € X:f(x) = I} and | f(x)\ < lforx E€ X¥ —S. 
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We say that such a function f peaks on the set S. 
The followings properties of A peak sets are elementary. 


Property 2.5—If S is an A peak set then S is compact and t (S) = S. 


Property 2.6—If S = A S, + ¢ where each S, is an A-peak set then S is also 
n=1 


an A-peak set. 
There exists f, € A such that 
= {x € X: f, (x) = 1} and | f, (x) | < 1 on 
X — S, for all n. 


ca Ln 
let f= 2 2"€ A. Clearly f peaks on S. 
n=) 
Definition 2.7—For f € A we define || f ||» = bad | f (x) | for any subset y of x. 


Property 2.8—Let S = a Sx 4 ¢ where each S, is an A- -peak set and Y bea 
closed subset of X, such that Sy ay Y = ¢. Then there exists f € A, such that || /|| = 
f|s = 1and||f Illy < 1. 

ProoF: Y. S=¢. HenceeYC S*= VU S;, where SS is open for all «. 

“cE A 


Therefore for some n 


Yas Y Sy implies ¥ A CO Sa ) = 9%. Now by 2.6 A Se is an A peak 
= = tm] 


set. Therefore exists f € A such that 


Aes , = RhE XS (x) = 1} and AG) Ls Lon Xess 


hse! im™1 


But 


SCO Se. 


{=I 


therefore f (x) = l on S. 


IfM=11¢@)1<1lonx— 4 Sa = 9 SS. 


i=1 im~1 i 
Since 


Ave ase 


i™1 ee 


lf (x)| <lony. Hence || f ||, < 1 
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by compactness of Y. 


The following theorem gives the relation between maximal A anti-symmetric sets 
and A peak sets. This theorem will be used to prove the main results of the paper. 


Theorem 2.9—Let K be a maximal A antisymmetric set. Then K U + (K) is the 
intersection of A peak sets that contain K U t (K). 


Proor : Let % be the family of all A peak sets that contain K U +(K). Since 
X is an A peak setand K U +(K) C X¥wehaveX¥€ & 
Thus G is non empty. 


Let Z= () S. Note that'Z = + (z), 
SER 


K U7t(K)C Z. We want to show that K U + (K) = Z. 
Claim 1\—f € Aand/f|z is real implies f| zis constant. (The proof of this 


claim follows closely the proof of a similar theorem for a complex algebra given in 
Burckel’]. 


If not, there exists f € A such that f | z is real and non constant. 


Since K is an A-antisymmetric set f | K is constant and hence f | + (K) = f | K is 
also constant. Thus f| K U +(K)isa constant saya. Since f| zis non-constant, 
there exists Z, € Z —(K U + (K)) such that f(z.) = b #a. Define go =1—( f—a)?. 


i -=.al, 
Then go € A. 
& (Z) C [0, 1] go | K U + (K) = 1 and g% (20) < 1. 
Let 
V, = {x € X: | 20 (x) | < 1 + 1/2"}. 
Then V,, is an open neighbourhood of Z = a) ss. 


SER 


Hence Ve CZ°= US. 


Thus {S*: S € H} is an open cover of the compact set Ve : 


Hence Vi C Si U S; U... US), for some mi.e., Ss S21... Sm CVn 
But SiN S, .. S,, = Sis an A peak set by 2.6 and K Ut (K) C S. 


Hence Z C S. 
Since S is an A-peak set there exists g, € A such that 


llgnll = 1 8, (S) = {1} 
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and 
IZnlle-y <1. 


By taking sufficiently high power of g, if necessary we may assume 


ed tae 
2” |lSnll * 
Thus we have |lg,|| = 1 g, (Z) = {1} 


IIgnlle-v < 
n 


and 
I wT} 
lSnllx-—Vn Fn an lio. : 
Let 
g=f = eee .-:(2) 
n= 2 
Then 
| g (x) | < | go(x) | for allx € X. wat) 


Ifx €  V, then | go (x) | < 1 + 1/2" for all n that is | go (x) | < 1. 
Thus by (3) 


lg(x)| <1 forallxE 2 V,, . (4) 
Now suppose x € (X — () V,). 
Since V, is a decreasing sequence (V,4 C V,) letting V) = X, there exists m > 1 
such that x € V,,_, and x & V,, then x & V, foralln > m. 
Thus ‘for = 11,22. c a —1 


| 80(X) 8 (OI S14 sey. (3) 


For n > m since x & V, 


| So (x) Bn (x) | < 1/2" < 1]/2™-2, ...(6) 
Combining (5) and (6) 


le@) 1< FF le) g(x)! 


m—1 


1 »S ] 1 
= SS l 
S (1 Qm-1 ) 2" Qm—1 Qn = | 
n=] 


that is 


lg(x)| <1 foralxE X-O)Y,. ...(7) 
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(4), (7) implies ||g\|,0 <1. 
Hence the function ery) € A peaks on the set S = g-! {1}. 


But go | K U + (K) = land g, | K U+(K) = 1 forall». Since K U t (K) C Zand 
&ilz= 1 


Hence g (x) = gy (x) 2 (1/2") &n (x) = 1 for all x © K Ur (K). 


Hence K U +(K) C S. 


This means that S is an A peak set containing K U+(K). Hence SE ZH. This 
implies that Z C S = g-} {1}. 


But 2) € Zand | g (zo) | < | 80 (zo) | < 1 by (3) which isa contradiction. 
This proves the claim. 


Now if f € A,f| zis purely imaginary implies f| z is constant and Z is an 
A antisymmetric set. But since K as well as + (K) are maximal A-antisymmetric sets 
(Kulkarni and Srinivasan®, Lemma 2.10) contained in Z we have K = Z = ¢ (RK). 
This proves the theorem. 


If not, there exists A € A such that h | z is purely imaginary and non constant. 
But then h? | z is real and must be constant by claim 1. We may assume A? lz=-l. 


Thus 4 (x) = +i for all x € Z. 


Let 
Z, = {x E€ Z: h(x) = 3} 
and 
Z. = {x € Z: h(x) = — i}. 
Then Z = 2; U Z,, Z, O Z, = ¢ and Z, and Z, are nonempty. + (Z,) = Z, 
and +(Z,) = Z,. Since h| K must beconstant,h| K =iorh|K = — i. Hence 


K & Z1 or K ce Zo. 

Claim 2—Z, is an A-antisymmetric set. Let f€ A and f| Zi be real. But then 
Jf | Z. =(f | Z,) is also real, and hence f | z is real and by Claim 1, f| z is constant. 
Hence f | Z, is constant. 


Now let f € A andf| Z, be purely imaginary. But then fh | Z, = if | Z, is real 
and hence constant. This proves f| Z, is constant. Thus claim 2 is proved. 


If KC Z, by maximality of K, K=2Z,,7(K)=Z,. If KC Z,,K = Z,, 
t (K) =Z,. 
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In either case, K U t(K) = Z, U Z, = Z. This completes the proof of the theorem. 


Corollary 2.10—Let A be a real function algebra on (X, +,)and K bea maximal 
A-antisymmetric set. Then A | K is uniformly closed in C (K) (see Theorem 3.6 
of Kulkarni and Srinivasan® for a different proof). 


Proor : For f € A ande > 0 wecan find F € A such that ||Fllx < Il /llx + €¢ 
and F = fon K. 


We may assume that || f||y < 1. 
Let 
V=(x € X:|f(x)| < Ifilx + 4. 


Clearly V is a neighbourhood of K U +(K). Hence by Theorem 2.9 and Property 2.8, 
there exists g € A such that |lg|| = 1 g | KU +(K) = 1 and |lgllx-v < 1. 


By taking sufficiently high power of g, if necessary, we many assume 
llglla-v < Ilfllx + «. 
Let F = fg € A. Clearly F = fon K. 
For 
x EV, | F(x) | = | f(x) Ile (*) | < Ilfilx + € as 
| f(x) | <Ilfilk +. 
Forx € X — V, 
| F(x) | = |f@) llg @) | < If llk + as 


tg (x) | < Ifilx + «: 
Thus 


|Fllx < Ilfllk + «. 
This proves the claim. 
Now let 
Ar ={fEA;f| K =0} 
and Q = A | Ax the quotient space. Q is complete in the quotient norm. 
For a typical element g + Ax of O, where g € A 
lig + Axll = inf lig + fl 
fEAK 


I 


inf {|All :h| K = g | K} 


= |lgllx by the above claim. 
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Thus A | K is isometric to Q. Hence A | K is also complete and hence closed in C(K). 


3. A CHARACTERIZATION OF C (X, +) 


We present a proof of the main theorem in this section using the properties 
discussed in Section 2. 


Theorem 3.1—Let A be a real function algebra on (X, 7). Let Re A = {Re/f: 
f€ A}. If Re A is uniformly closed then A = C (X, 7). 


We shall need the following Lemmas in the proof of the above theorem. 


Lemma 3.2—Let A bea real function algebra on (X, 7). If Re A is closed and 
A is antisymmetric then X is a singleton set. 


Proor : Since A is an anti-symmetric algebra, by Definition 2.3 we have the 
following : 


(i) f € A and fis real implies f is constant and 
(ii) f € A and fis purely imaginary implies that f is constant. 


Let x) € X. Since A is antisymmetric, for every u € Re A there exists unique 
f € Asuch that u = Ref. ' 


Define T (u) = f. Then T: Re A => 4 isa linear map. Let u, € Re 4,u, > u 
and T (u,) = f, >f. Then clearly |Re f, — Ref lloo <Ilf, — fll > 0 


that is uv, ~ Ref and hence u = Ref. 


Thus T(u) =f. This shows that T has a closed graph. Hence T is a bounded 
operator. Suppose X contains more than one point. 


Let g =u + iv € A bea non-constant function such that g (xo) = 9. 


Let 
i= ae iaied ee st 
— |p|] < » S Ill 
Since g is non constant V # 0. Hence there exists x, € X such that | V(x) | = IVI. 


Let Z, = g(x:). Then Z, and Z, are boundary points of R. 
Let 

N > Oand Ry = a Sees it 

and—N<y€Nn : 

[Note that R and Fw are symmetric w.r.t. real axis, that is w € R (respectively Rw) 
if and only if # € R (respectively Rw]. There exists a conformal map ¢: R° ~Ry 
such that 

¢ (z;) = iN and ¢ (z;) = — iN. 
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This map extends to a homeomorphism of boundaries of R and Ry. 

Thus ¢ is a uniform limit of polynomials in R. 
Define |: R > Rn as ¥ (z) = $ (z) + $(). 
Then ¢ is a uniform limit of polynomials with real coefficients because if {p, (z)} isa 
Sequence of polynomials converging uniformly to ¢ (z) on R then gq, (z) =p, (z) 
+ pa (z)isa sequence of polynomials with real coefficients converging uniformly to 
(2). 

Therefore h = ) (g) € A by spectral mapping theorem‘. Now T (Re A) = h, 


h(x) = $ (gx) = ¥(z) = 4) +4) 
= iN+ (—iN) = 2iN. 
Thus we have |/Re Al] < |lul| 


[hl] > N. 


Since N is arbitrary this is a contradiction to the fact that Tis a bounded linear 
map. Therefore X is a singleton set. 


Lemma 3.3—Let A be a real function algebra on (X, +t) and suppose Re 4 is 
uniformly closed. Let S be a maximal A antisymmetric set. Then 


Re Aj’ {Re f|s:f € A} is also uniformly closed in C (S). 
Proor : Claim 1—If f€ A ande > 0 wecan find F © A such that [lRe Fllz 
< [Re f lls + 2€ 
and Re F = RefonS. 
We may assume that || /|| < 1. 
Let 
= Ww=s+itEe C:|w| <1 and|t| <6. 
Let 
U={zE€C:|z| <1}. 


(Note that both 0 and U are symmetric w.r.t. real axis) and v be a conformal map of 
U on © with y (0) = 0 and v(1) = 1. 


Now v is analytic in U, 


Define »: U+> Qas 


p (z) = (v (z) mot) 
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Then up is analytic in U, u (0) = 0 and p (1) = 1. 

Choose 6 > 0 such that | z | < 8 implies | » (z)| <. 

Since S is a maximal A antisymmetric set, by Theorem 2.9, S U +(S)= © Aw 
where each Aq is an A-peak set. eS 
Let 

V={xE X: | Ref(x)| < |[Reflls + ¢. 
Then V is a neighbourhood of S. 


Since Re f (x)=Re f (t (x)) for every x € X, we have V=+ (V). Thus S U + (S) C V. 
Therefore by 2.8 there exists g € A such that ||g|| = 1 


g|S U7 (S) = 1 and |lgllx_v < 1. 
Choose a positive integer n large enough such that 
llg"Ilx-u < 8, 


Let h = »p (g"). Thenh € A by the spectral mapping theorem for real Banach 
algebras‘. 


h=1onS U t(S) 

|Imh| < € on all of X. 

Also | Reh | < ¢€on X¥ — Uand|ReA|<lonallof X¥. LetF=fhe A. 
Re F = Ref Reh — ImfIm A. 

h = 1 on S implies Reh = 1, Imh=OonS 
Re F = Ref onS. 


On 
U| Re F| <|Ref||ReA|+|Imf||ImA | 
< (IRefils +). 14+ « 
On 
¥—U,| ReF| a «+ ¢, 
Thus 
[Re F |lx < |IReflls + 2¢ 
and 


Re F = Ref on S. 


To prove Re As is closed : Let Rs = {u € Re A:u(S) = 0} ReA with the 
supremum norm is a Banach space and Rs is a closed subspace of Re A. 
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Consider the quotient space 
QO = Re A/Rs = {Ref + Rs: f € A}. 
Then Q is complete in the quotient norm 
[Re f + Rs|| = inf {[Re Fl]: F € A and Re F = Ref on S} 
= ||Re f |ls by claim 1. 


Thus Re As is isometrically isomorphic to Q and hence complete in the maximum 
norm on S. 


This proves that Re As is closed. 
We make use of the following theorem in proving Theorem 3.1. 


Bishop’s Theorem’—Let A bea real function algebra on (X, t). Suppose f € C 
(X, t) is such that f | x € A | .for every maximal A-antisymmetric subset K of X, then 
FEA. 


Proof of Theorem 3.1—Let K be a maximal A-antisymmetric set. We shall show 
that K is singleton. 


Denote the restriction algebra 4 | x by Ax. Note that Ag is uniformly closed 
in C (K) by Corollary 2.10. By Corollary 2.11 of Kulkarni and Srinivasan®, either 
K = <'(K) or. K()-+ (KV = ¢. 


Case (i) K = ~(K): In this case, Ag can be tegarded as a real function algebra 
on (K, + | x). Obviously Ax is an antisymmetric algebra. Re Ax is uniformly closed 
by Lemma 3.3. Hence by Lemma 3.2, K is singletone. 


Case (ii) K (1. t (K) = ¢: In this case it can easily be shown that there exists 
h € A such that h| x =i (see Theorem 2.17 of Kulkarni*). This meansi € Ag. 
Thus Ax can be regarded as a complex function algebra on K. As before, Ax is an 


antisymmetric algebra and Re Ax is uniformly closed. Hence K is Singleton by the 
main lemma in Hoffman and Wermer?. \ 


Thus every maximal A-antisymmetric set is Singleton. Hence by the above 
analogue of Bishop’s theorem A = C (X, =). 


Corollary 3.4—If Re A D Cr (X, t), then A = C (YX, vt). 


Corollary 3.5 (Stone-Weierstrass)—Let A be a real function algebra on (X, 7) 
such that f € A for allf € Athen A =C (Xo): 


> 


Proor : Since A is closed under complex conjugation Re A is the collection of 
all real valued functions in A. Hence Re 4 is closed. 


Remark 3.6 : Mehta® has defined a real function algebra on a compact Hausdorff 
space Yas a uniformly closed real subalgebra of C (X) that contains the constant 
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function | and separates the points of ¥. If A is such an algebra and Re A is closed, 
then what can be said about A? The following theorem answers this question. 


Theorem 3.7—Let X be a compact Hausdorff space and Aa uniformly closed 
real subalgebra of C (X), which contains (real) constants and which separates points of 


X. If Re A is uniformly closed in C (X), then there exists a closed subset Z of XY such 
that 


A= {fe C(X):f| zis real}. 
Proor : Let Z = {x € X: f(x) is real for all f€ A}. LetY =X x {0, 1}. 


For every z € Z we identify the points (z, 0) and (z, 1) and denote the space obtained 


by this identification by XY. (We may consider Y as the space obtained by pasting 
together two copies of X at the points of Z). Define +: Y > Yas + ((x, 0)) 


= (x, 1) and « ((x, 1)) = (x, 0) for every x € X. Let be the map induced 
by t on X. It is easy to see that X is a compact Hausdorff space and + is an involu- 
toric homeomorphism on y. We may regard X as a subset of Xx . Then every point 
of Z is a fixed point of B 


Now Let f € C (X) and f be real on Z. 
Define 


f on X as f(x, 0)) = f(x) 


F(x, )) =F) 
for allx € X. 


Since for x € Z, (x, 0) and (x, 1) are identified and f (x) real, f is well-defined. 
Clearly f € C(X, t) and f| X¥ =f. 


Let pue WG] € A}. Then it is routine to check that A is areal function algebra 


on (x, «) and Re A is uniformly closed. Hence by Theorem 3.1, A = C(X, 1). This 
in turn, implies that A = {f € C (X): fis real on Z}. 


Remark 3.8 : As special cases of the above theorem, we get : 


If Z = ¢, then A = C(X) andif Z = X, then A = Cr(X). In particular, if 
for every x € X there exists f € A such that Im/f(x) # 0, then Z = ¢ and hence 
A = C(X). This is a stronger version of the classical Hoffman-Wermer theorem. 
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co 
Let f(s) = = an eXp (s\n) represent an analytic function in the half plane 
n= 


Res <a, (—2 <a <oo), Let Ey (f, 8) be the error in approximating the 
function f(s) by exponential polynomials in the half plane Res <8 < « In 
this paper, the lower order and lower type of f(s) have been characterized 
in terms of the rate of decrease of E, (f, 8). Further, the order, type, lower 
order and lower type of f(s) have been obtained in terms of the ratio of con- 
secutive error terms En; (f, 8)/En Cf. 2). 


§ 1. Let 


f(s) = z a, exp (SA,) SHER) 


where 0 < A, < A,,; > 00, 5 = ¢ + it (o and ¢ being real variables), [a,];° is a sequ- 
ence of complex numbers, be a Dirichlet series. We assume that 
lim inf(A,41 — A,) = 90> 0. (Ua) 
n~>o 
Then we have 
lim sup (n/A,) = D < o9. ples) 


n->o 


If the series (1.1) converges absolutely in the half plane o = Res <« then */ (s) 
is an analytic function for o < « and 


ay 
« = lim sup slog Ss ; ...(1.4) 
n-> oo n 
= |. u. b. o + it) | . The growth of the analytic 
For any > < 4, we set M (s) ar ata 


function f (s) is studied through the parameters such as order, type etc. Thus /f (s) is 


said io be ? of order P if 


; loglog M()_—ss ep OS p<. ..(1.5) 
lim sup fog [i — exp (« — 2)] : 


oa 
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Further, if 0 < p < co then the type 7 of f(s) is defined as 


lim sup OR ae (2 


poe nied Deis AL PR gy OVE age y BP ae tho) 
oa [1 ae exp(o os a)] , 


Likewise, the lower order A and lower type (0 < P < oo) of f(s) are defined as*’5, 


log log M (c) 





lim inf ——>—">—— —""___ ~) 0<A<p<c eA teal 
—log [1—exp («—<]) 
lim inf STE A elas EO eee ce (28) 
ona) Ll ee OxpASiees a)? 
Let »(s) = max { | a, | exp (oA,)} be the maximum term inthe expansion (1.1) and 
nF 1 
An(e) = max {u (c) = | a, | exp (cA,)} denote the rank of the maximum term. For an 


analytic function f (s) of finite order, it is known Nandan® that log M (c) ~log # (cs) as 
o — a. Hence M (s) can be replaced by p (c) in the definitions (1.5) to (1.8). Nandan? 
proved that if f(s) is of lower order ? then 





a log Ane) 

1 + A = lim inf — : --.(1.9) 
o> lor [\— exp (¢ — a)] 

However, it was shown® by means of a counter example that the result (1.9) does not 


always hold. We shall call an analytic function f(s) to be of “admissible class” if 
(1.9) holds. 


Let D. denote the class of all functions given by (1.1) which are analytic for 
¢ < «and let Ds, — co < B < « < co, denote the class of all functions given by (1.1) 
which are analytic for o < 8, The notion of approximation of an analytic function by 
means of exponential polynomials was introduced by Nautiyal and Shukla’. In the 
series (1.1), ifa, = Oforn > k and a, <0 then f(s) is saidto be an exponential 
polynomial of degree k. The class of all exponential polynomials of degree almost k is 
denoted by 7%. Following Nautiyal and Shukla’ we define 


E, (f, 8) = inf lf — plie,n = 0, 1, 2,... ...(1.10) 


PET y 


where || f — pile = 1. y b. | f(B + it) — p(B + it) |. Nautiyal and Shukla obtained 
—oo x. f< 


oo 


the characterizations of P and Tin terms of the rate of decay of the approximation error 
E, (f, 8). Hence (Nautiyal and Shukla’, Theorem 2 and 3) we have 


Theorem A—Let f(s) € Dz be of order P and — co < B< «4 < ©. 
Then 


ro. log* log* [E, (f. B) exp {(« — 8) Anat] 
= | 
ET. pees 2 log heeeceloee eset (7, B) exp {(@ — 6) Anal] 


A018 


APPROXIMATION OF ANALYTIC FUNCTIONS 169 


Theorem B—Let f(s) € Dx. and —ce < B <a< oo. Then f(s) is of order 
P (0 < P < co) and type 7, if and only if 


P+1 T 
(2 ) (PT) = “lim ‘a [logt{E, (f, den — B) A,.1)}]et? 


| (142) 
when the right hand side of (1.12) is a non-zero finite number. 


In the present paper, we obtain the characterizations for the lower order A and 
lower type ¢ in terms of the decay of the approximatione rror E, ( f, B). We also obtain 
these characterization in terms of the ratio of consecutive errors E,-, (f, 8)/E, (J, 8). 
To avoid repititions, we shall assume throughout that f(s) € D. and the condition 
(1.2) is satisfied. Further, we shall write E, for the error term E, ( f, 8). 


§2. We now prove 


Lemma—Let fe (s) = = {E,-1 exp (—BA,)} exp (sA,). Then fe (s) € Dz. 


Further, 
T = T (fe), t= t (fe), 0<P<co 


Proor : That fa (s) € Ds follows from (Nautiyal and Shukla’, Theorem 1] and 
(1.4). Consequently, we have Pp = P ( fe) from (1.11) and Theorem | of Juneja and 
Nandan!. From Lemma 1 of Nautiyal and Shukla’, we have 


E,-1 < K M (s)/exp ((¢ — B) A,), 2 = 1, 2, «.. ...(2.2) 


(0) 


where K is a constant independent of n and «. Hence we have 
(9, fa) < K M(s). 
Also from (3.6) of Nautiyal and Shukla’, we have 
M (cs) <| a | + 2M (c, fp). ...(2.4) 


Hence from (2.4) we have \ < (fg). Further for 0 < P < oe, (2.3), (2.4) and the 
definitions of T (fe), T etc. lead to T= T (fa), t = t ( fa). This proves the lemma. 


In the following results, we denote 
logtx = max (0, log x). 
Theorem 1—Let f (s) be of order P > 0 and lower order ,0<A<o. 


Then 


eee One ee, (2.5) 
log Angi — logt log* [E, exp {(@ — 8) Ans] 





1+A > lim inf 


n->e 
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Further, if Y (n) = log (Ey-1/En)/(An¢i — An) forms a non decreasing function of 1 for 
all large values of n and fp (s) belongs to ‘admissible class’, then 


bats, ¢ log A, 
Eh Le ented Nad En Es ...(2.6 
1 +A lim ‘inf log An+1 — log* log™ [E, exp {(« — B) Anti] ee 


n—> oo 


Proor : Let the right hand side expression in (2.5) be denoted by 9. We assume 
that 1 < @ < oo, Fora given e > 0 such that 0 < « < @, there exists integer N=WN (e) 


such that 
log Ey + (% — B) Ansa > Ants Ani 1O-9 n > N (6). 
From (2.2), we have for all « close to «, 
log M (c) > log E, + (6 — B) Ak+i + O(1), & = 1, 2, ... ESS 
Hence we have 
log M (¢) > Angi (An) 2MO-9 + (6 — «) Ani, 2 > N (e). 
We choose a sequence [s,] by defining | 
a — G, = (Ie) (A,)- 1-9), n > N. 
It is evident that 6, < « ando, > « asn + oo, Foro, € « < a,4;, we thus have 
log M (6) > Anya (An) OO + (6, — &) Ang 
=(1 = e7*)Angr (An) O79 
= e(1 —e4)(a — o,) [e (@ — o,,,)] “9 
> (e — 1 ) (a — a)!-O-9 e-(o-«). 
Therefore 
log log M (s) > O(1) + (1 — @ + €) log (a — a). 


It can be easily seen that 1 — exp (« — «) & « — caso > «. Hence on proceding to 
limits, we obtain 


Ws log log M (cs) : 
A=] = 
wes Gone reg exp —(¢ — a)] 2 @— I, 


The result is trivial if@< 1. If @ = oo, then the result follows on taking an arbitra- 
rily large number in place of (@ — e) and Proceeding as above. This proves (2.5) 


To prove (2.6), we consider fj nese 
( ider fp (s) 2 b, exp (sA,,) where b, =E,-, exp (—BA,). 


Then 
— 1081 Olbnss, log (Ey-1/E,) 
§ () Ans Any Or Lee = iB ath (mn) Be were et 
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Hence ¢ (n) is a non-decreasing function of m whenever / (n) satisfies the condition. 
Applying Lemma 3 of Nandan‘ to the function fe (s) (under the assumption that (1.9) 
is satisfied for fp (s)) we obtain on using (2.1), 


on log A,,-1 
1] at ey ee = ee Cd ee eye re 
+ Oye eee UE toe a Jap oat [Esc sip Tec BL 


This proves (2.6) and the proof of Theorem | is complote. Next we prove 


Theorem 2—Let f(s) be of order P, 0 < p < co, and lower type t,0 < t € o0, 
Then 





\erl 
i ae iF (pt) > lim inf A, 


p N->co 


| ee {E, exp A — B) Ang) ke 
Ani 


-..(2.9) 


Further, if | (nm) as defined before, forms a non-decreasing function of for all large 
values of n, then 


An+1 


“d — +1 
(f+) \P" (et) < tim inf Aya. pews B) Anas} et. 


noo 


...(2.10) 


Proor : Let the right hand side expression of (2.9) be denoted by 4, 0 < A<co. 
Then for a given e, arbitrarily small and satisfying 0 < « < A, we have 


+ 
A Hee Ex EC) ea] SAH en — N= N:(4) 
n+1 


or 
log E, > Anyi [{((A — OfAgPMPTY — (2 — B)),n > N. 


On using (2.7) we have 
log M (c) > Aner [{(A — €)/A,} et) — (a — )] + O(1),n > N. 


We define a sequence [c,] by putting 


oe 
a —s6, = ( an 7) (At), or =e NV + 15.0.: 


Then for sc, <o < 6,4:, we have 


log M (0) > Anyi [(A — €)fAq]CP*9 — (a — Gy) Ang, 2 > N 


it real A Vi a —o)-*, as in Theorem 1. 
= ot sai ) (A — «) (a —0)-*, asi 


Hence we get 


Ber. log M (c) Z( ee a 
t = lim inf (= expe — a)}? 2 p p+ il 


a--a 
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and (2.9) follows. The proof of (2.9) when A = oe can be completed as before. 


On using the condition on the function ¥ (n), we obtain for analytic function 
fe (s). [Nandan®, Nandan and Srivastava® (Th.*)] 
log* {E, exp ((% — ced 


nm Anta 


(° + ] ie pt (fp) < lim inf 


Inequolity (2.10) now follows on using (2.1). This proves Theorem 2. 
Combining (2.9) and (2.10), we obtain. 


Corollary—Let f (s) € Dx be of order p, 0 < Pp < ce, and lower type 4,0 <t 
log (E,—/E,) 
A 


forms a non decreasing fun- 
Anti— n 


<S oo. If A, ~ Angi asn > co and (nm) = 
ction of n for all large values of n, then 


(se i (rt) = limint | 208 CEs eee (BAD) | 
p n 


Neco An+i 


a VEL E 


§3. In this section we obtain the characterizations of P, A, 7, and tin terms of the 
ratio of error terms (E,,/E,). We prove 


Theorem 3—Let f (s) € Da be of order p and — co < B<a < oo, 


If {(n) = We Cells) forms a non-decreasing function of n for all large 
n+ n 


values of n then 


oa logAn41 
Pt, sim SUD 16g (,a-A.) log [=A (4, aang (Ege eal! 


...(3.1) 


Proor : Let the expression on the right hand side of (3.1) be denoted by B+ 1. 
We first assume that B < oo, Then fore > 0 and arbitrarily small, 


- 


log Ane1 
log (nea — An) —log[(@—8) (nui — An) + log (EE) 2 TET 
n> N¢e) 





or, 
log (En/En—1) < (Ansa — An) Angi) 8449 — (a— B) ngs — An) 2 > N (0. 


Substituting nm = N + 1, N + 2, ..., kin the above inequality and adding, we obtain 
k k 
log (Ex]Ew) <_& (Anti —An) Angi) 49 — (a—B) E (Ang, — Aj). 
n=N+1 n=Nt+1 


Putting c = 1B + 1+ and rearranging the terms on the right hand side, we obtain 
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k 


log Ex < AL — Ana ANS. — > A, (ax, — a;*) 
n=N+2 


+ (a — B) (Aw4i: — Akai) + log Ew. 
Let us put 7 (¢) = X, for A, for A, St < Angi and F(t) = 7, 


Then 
Mets 
lorsep ‘log’ By <P AY” Ay, Ays = | n(t) dF (8) + (a—B) Aner —Aeaa) 
AN. 
Aa 
Sine yeaa c| Boe f-odt =" e BAGS. 
“N41 


Since n (t)/t < 1 we obtain 


flog Ex + (@ -) Mul <O) +ALL + TOGA RDN 


1 -c 
=O (1) a ee 


Precceeding to limits as kK — ©, we obtain 


log* log* [Ex {exp ( — B) Ax,}] — BB 


lim sup logue Wei 


ko 


and in view of (1.11), we obtain p/P+-1 < B/B + 1. Hence we have p < B. 


The above proof obviously holds for B > 0 only. However, if we assume B < 0 
then for suitableec > 0,1 —e <Q. 


Then 
log* {Ex exp ((% — B) Ax4i1)} < O (1) 
and hence from (1.11) it follows that P = 0 while we have assumed that 0 < Pp < oo. 
The inequality obviously holds if B = °°. 


To prove the reverse inequality, we use the condition given for J (n). Then for 
all n > no (fixed), we can write 


log (En [Ex) = = (s) Are — ds) 


< p (n) (n4i my + Fa 
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or, 
log* E, 
(a — p) + EEE) (a — B) + Reta O (1), n > mo. 
Pat eas Gel 

Hence 

li log An+l 

nrc | TOR Ans: — An) — log [(% — B) (nga — An) + log (Eql Ey) 

< lim sup 10g Ant 





noo log Anis — logt [logt E, + (% — B)Angi) * 


Hence B + 1 < P + | in view of (1.11). This proves that B = P and hence Theorem 
3 follows. 


Theorem 4—Let f(s) € Dx be of order Pp > 0 and lower order \, 0 < A < ©, 
Then 


hor loga, 
Poe Met ant Toet Wasa @o By Oca eh ae ee 
...(3.3) 


Further, if ) (n) = 108 (BacsiE) forms a non-decreasing function of n for all large 
n+1 — 4n 


values of n and fp (s) belonge to the “‘admissible class” then equality holds in (3.3). 


Proor : Let us put the right hand side of (3.3) equal to 6. We first consider the 
case when | < §@ < oo, For arbitrarye, O<e < 6, we have 


log (En/En-1) > (Ansa —An) {Ay “I — (a — BY}, n > N= NCC) 


Proceeding as before we obtain 


- /(€@) 
log EE. = log En + > (Ak41 = An) fA,’ —(a — B)) 
k=N+1 : 


> log Ew + [A7@-o) — (q — B)] (A 


Rely. An)s 


since A, is an increasing sequence, Proceeding an in the proof of Theorem 1, we can 
prove what 1 + A> 4. This inequality obviously holds if 0 Se at 


For @ = co, 
above method gives (3.3) if we replace (9 — e) by any arbitraily large num 


ber. 


To complete the proof, we use the condition on } (n). Then from (2.8) and 
Corollary of Nandan‘ (p. 1367), we get for the analytic function fe (s) satisfying 1.9), 


1 9 (fo) = lim inf oD oe 
fe) = Troe oa? Gaps) Sea (n+/bn)i] 
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where b, = E,-1 exp (— BA,). In view of (2.1) and (3.3) we now get | + A = 0. This 
proves Theorem 4. 


Lastly we prove 


Theorem 5—Let f (s) € Dx be of order P, 0 < p < ©, type T and lower type ¢, 
O<t<T< oc, 


(i) We have 


= 


PT Ke = are. Io (PT) ...(3.4) 


where the right hand inequality of (3.4) holds if & (n) forms a non-decreasing function 
of n for all large n and 


6 = lim sup An4, 
n->oo 


log (E,/En-1) + (% — 8) (Angi — An) tis 
[ Ant Aa TAg 


(ii) If A, A,,, asin > cc, then 


“i — — Atl. 
pr > lim inf A, [ a DOES led Jd Wes on e oe 


egy nis ss 


ll co 


Proor : The right-hand inequality of (3.4) follows immediately on using (3.2) 
and (1.12) of Theorem B. The proof of left-hand inequality of (3.4) i.e. pT < @ 
follows on the lines of proof of first half of Theorem 3. Hence we omit the proof. 


To prove (3.5), we denote the right-hand side expression of (3.5) by g. Then 
for arbitrary «, 0 << « <q, we have 


log (E,/En—1) + (@ — B) (Anza — An) > (Ants — An) ( i= jue 
n> N(e). 


Now 


b 
log E, = ne log (E,/E,-1) — log En 


> > Gen ay (SY - 5 (« — 8) (nti ~ Aa). 


n=N+1 n=N+1 
7 log En. 
Putting c = 1/(P + 1) and rearranging the trems on right-hand side we get 


k 
log Ex > (q — ©)° Det Ay — Ayes >. An ( a te as*)] 


n=N+2 


— (a — B) Anyi + O(1), K > N. 
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We put n (t) = A, forA,-1 << t SA,, F(t) = t-%. Then 





ay 
log Ex > (q — €)° (Asi Ag” — Ana — | n (t) dF (t)] — (@ — B) Ans 
ANG? 
at) (a) 
or 
Ay 
log Ek + (2 — B) Ain > O(1) + @— OF Bann dy? +e "OF ay 
Aves 
Since n (t)/t > 1, we have 
log Ex + (4 — 8) Acar > O (1) + ( — © [eer dy + 5 
(q—«r Ay 


(A, °— Av) > OM) + Te fe 


since A, is an increasing sequence. Hence on proceeding to limits, we obtain 


lim inf af w08" Fs + (@ — 3) Aen IF . (a) q. ...(3.6) 


Since A, ~ A,,1 as n —> ec, combining (3.6) and (2.9) we get (3.5). This completes the 
proof of Theorem 5. 


—- 
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In this paper a characterization for the spatial numerical range of a normal 
operator on a smooth reflexive Banach space, to be closed and convex, is 
given (Theorem 1). This generalizes the Theorem of Ching-Hua-Meng’, for 
normal operators on Hilbert spaces. A few more results concerning the 
spatial numerical ranges of convexoid and isoabelian operators are also 
obtained. 


Let X be a complex Banach space and B(X) the Banach algebra of bounded 
linear operators on X. For T € B(X), the spatial numerical range V (T) and the 
numerical range W (7) of T are given by 


Vy (f(t) xe Xf e'xX*, Hl St =F (x) 
and 

W (T) = {{Tx, x] : |Ixll = 1} 
where [.,.] is a consistent semi-inner-product (s.i.p) of Lumer®. The sets V (T ) and 
W (T) are neither closed nor convex in general. It is easy to see that V (T) is the union 


of all the numerical ranges W(T) corresponding to all consistent s.i.p’s on X. 
However, V (T) = W (T) in case of smooth space. 


T € B(X) is said to be hermitian if W (T) is real or equivalently V (7) is real. 
If 7 = H+ iK, where Hand Kare hermitian with HK = KH, then T is calleda 
normal operator. The numerical range of T w.r.t the Banach algebra B(X) is 
defined by 
V (B(X),T) = {F(T): FE B(X)*, IF ll = F (I) = I}. 
It is easy to see that W(T) C V (T) C V(B(X), T). 


An extreme point of a convex set S C C isa point p of S which does not belong 
to any line segment joining two points of S which are different from p. 


The sets « (T), Co o(T), % (T), % (T), or (T), 54 (T) and T (T) respectively denote 
the spectrum, convex hull of the spectrum, continuous spectrum, point spectrum, residual 
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spectrum, approximate point spectrum and the compression spectrum of 7. An operator 


T is called convexoid if V (JT) = Co«(T). A normal operator on a Banach space X is 
convexoid. Indeed, if J is a normal operator, then V (B (X), T)=Co o (T) (Bonsall'and 


Duncan’, Th. 14, p. 54). Since Cos(T) C V(T) for any T € B(X) (Bonsall and 
Duncan’, Th. 4, p. 22) and V (T) C V (B(X), T), we have V (T) = Coc (T). 


Throughout this this E denotes the set of allextreme points of V(T) of a 
convexoid operator 7. 

Definition— A Banach space X is said to be smooth if for each unit vector x in X, 
the set D(x) = {f € X*: || fll = f (x) = 1} contains only one point. 

Theorem |—Let T be a normal operator on a smooth reflexive Banach space XY. 
Then V (T) is closed and convex if and only if EM 9. (T) = ¢. 


ProoF : Since T is normal, V (T) = Coc (T), and therefore E C o (T). Suppose 
E (\ o¢(T) = ¢. Then as o, (JT) = ¢ (Mattila®, Th. 4.7), E C op (T). Since V (T) is a 
compact convex subset of C, we observe that V(T) = CoE. As Coo,(T) C V(T) 


(Bonsall and Duncan’, Th. 3, p. 21), it follows that V(T) C V(T). Thus V(T) is 
closed and convex. 


Conversely, suppose V (7) is closed and convex. Since V(T) C V (°C) 
(Bonsall and Duncan’, Corollary 3, p. 12), V (7) = V (T*). T* is normal because T is 
normal. As X is smooth and reflexive, X* is strictly convex, and therefore 


E = EQ V(T) = E(\ V(T*) C oy (T*) (Mattila®, Th. 7.2). 


Now the reflexivity of X and X* impliles that o, (T) = o, (T*) (Goldberg’, p. 71) and 
or (T) = ¢ = a, (T*) (Mattila®, Th. 4.7). Therefore, 6, (T) = o, (T*), and hence 
ECo,(T). Thus Eo. (T) = ¢. 


We note that Theorem | is an extension of the Theorem proved by Ching-Hua- 
Meng? for normal operators on a Hilbert space. 


- 


Proposition 2—Let T be a normal operator on a separable reflexive Banach 
space X. If E (\ o, (T) = ¢, then E is countable. 


ProorF : By the first part of Theorem 1,ECoa,(T). Since X is separable, o,(7) 
E is countable by (Mattila®, Cor. 3.10). 


Proposition 3—Let T be a hermitian operator ona reflexive Banach space. If 
E (1) « (T) = ¢, then there exists an eigenvalve « such that | «| = ||TI}. 


PROOF Since V (T) = Co a (T), ECo(T). As T is normaloid (Sinclair'*, Prop. 2) 
there is an «in o (T) such that [« |] = ||7]], and hence « is in E. By the first part of 
Theorem 1, E C a, (T). 
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Theorem 4—Let T be a convexoid operator on a Banach space. Then 
(ii) If E C o,(T), V (T) is closed and convex. 


. PROOF : (i) Since V(T) = Coc(T), EC o(T). Aso(T) C V(T), Ecanot be 
in the interior of 5 (T). Hence E C 80(T) C on (T). 


(ii) Since Cos(T) = V(T) = CoE and Coc,(T) CV (T) (Bonsall and 
Duncan’, Th. 3), it follows that V (T) is closed and convex. 


Remark 1; V(T) ~ V(T*) in general. However they are equal for a convexoid 
operator J on a Banach space with E C «,(T). This follows from (ii) of Theorem 4 
and Corollary 3 of Bonsall and Duncan? (p. 12). 


Remark 2: It would be interesting to know whether the converse of (ii) in 
Theorem 4 is true. However, we prove 


Theorem 5—Let T be a normal operator on a strictly convex Banach space. Then 
a necessary and sufficient condition that V (T) be closed and convex is that E C 4p (T). 


Proor : Since a normal operator is convexoid, the sufficiency follows from (ii) of 
Theorem 4. Also the condition is necessary for, if V (T) is closed, then 


E = EQ V(T)C a, (T) by (Mattila’, Th. 7.2). 


A non zero vector x in a normed linear space X is said to be orthogonal to y in 
X w.r.t a compatible s.i.p. [.,.] on XY if [y, x] = 0. This orthogonality is not symmetric 
in general. 


Lemma 6—If M is a proper closed subspace of a reflexive Banach space X, then 
there is a non-zero vector y in ¥ and a compatible s.i.p[.,.] on X such that y is ortho- 
gonal to M. 


This is our Corollary 9 of Puttamadaiah and Huchegowda”’. 


Theorem 7—If X is a reflexive Banach space and T€ B(X ), then ' (T) C V (T); 
in particular o, (T) C V Ch) 

Proor: Let » € I(T). Then (7 — p/) X is not dense in ¥. Then by Lemma 6, 
there is a compatible s i.p [.,.]on X anda non zero vector z in Y such that [((T — p/) 
x, z] = Oforall xin X. In particular, for ¥ = z,h = [Tz, z] with ||z|| = 1 and hence 
p € W(T) CV (7). 

Mattila’ (Th. 7.3) has proved that if T is a normal operator on a smooth reflexive 
Banach space, then EMV (T) C % (T); consequently E — V(T) C % (T) because 
o,(T) = ¢ (Mattila®, Th. 4.7). We now extend the latter result for convexoid operators 
on a reflexive Banach space. 
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Theorem 8—If T is a convexoid operator on a reflexive Banach space X, then 
E—V(TM)C «6 Wa 


Proor: Since V(T) = Coa(T),E Ca(T). As o,(T) C V(T) always and 
o, (T) C V (T) by Theorem 7, it follows that 


E— VG) oC) F(T) Cog l7), 


For each x in a normed linear space X, there is an x* in X¥* such that x* (x) 
= ||x||? and ||x*|| = ||x|]. Let ¢ associate each x in X to exactly one such x* in X* and 
ax to ux*. ¢ is called a support mapping. There are infinite number of such 
mappings unless the space is smooth. Every support mapping ¢ defines a compatible 
s.i.p [.,.] on X if we set ¢ (x) (vy) = [y, x]. An invertible operator T on a Banach 
space X is said to be is so-abelian if there is asupport mapping ¢ such that ¢7=T*"'¢, 
or equivalently (Tx)* = T*-! x* or [Tx, »] = [x, T-! y] forall x,y in X where [.,.] 
is a consistent s.i.p on X defined by ¢ (x) (y) = [y, x]. 


We note that the following three conditions are equivalent for a bounded linear 
operator T on a Banach space. 


(1) T is iso-abelian 

(2) T is an invertible isometry 

(3) T is invertible and ||T|| = ||7-|| = 1. 

(1) <> (2) (Koehler and Rosenthal‘, Cor. 1). It is easy to see the equivalence of 
(2) and (3). 

Theorem 9—If X is a reflexive Banach space and T € 8 (X) with 0 & V (7), then 

(i) Tis1 — 1 
(ii) TX =X 
(iii) If T is an isometry on X, then T is iso-abelian. 


PRooF : (i) Suppose x is a unit vector such that Tx = 0. Then by Hahn Banach 
theorem there is an f € X* such that fll = 1 and f(x) = 1. Now0 = f (Tx) 
€ V(T), a contradiction. We note that X need not be reflexive in this case. 


(ii) Suppose 7X + X¥. Then by Lemma 6, there is a compatible S.i.p [.,.] on 
X anda non zero vector y such that (Tx, y] = O forall xinX. In particular 


(Ty, y] = 0. Since we can take |ly|| = : 
: yl—=10EW(T)CV ee 
proves (ii). (T) C V(T), this contradiction 


(iii) Since T is an isometry, TX is closed ss : 
. ee , , and hence by (ii) 7 is ont it i 
an invertible isometry. Thus T is iso-abelian. Cy pate ang eee 
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A point » € a0 (T) is said to be a proper boundary point of o (7) if there exists 
abounded sequence (u,) in P (7) such that (Wn — #) (4, — T)-| > 1 asn > 9, 
where P (T) is the resolvent set of T. The set of all proper boundary points of o (T) 
is denoted by p, (7). 


We observe that if T is iso-abelian, then so is 7*. 


Theorem 10—If T is an iso-abelian operator ona reflexive Banach space X, then 
o,(T) = ¢. 

Proor : If » € o(T), then |u| =|I7|| and so pe OV (B(X), T). Since 
o(T) (1 BV (B(X),T) C p, (T) (Mattila’®, Lemma 1), » € P, (T) or equivalently 
0 € P,(T — pl). Since X is reflexive, ¥ = Ker (T — pl) @ (T — pNDX by (Mattila®, 
Cor. 4.5). This shows that o, (iy 1o. 


Corollary 11—If T is an iso-abelian operator on a reflexive Banach space Y, then 
op (I) = a, (T*). 

Proor : Since 7* is iso-abelian, the reflexivity of X and X* implies that 
o¢(T) = , (T*) (Goldberg’, p. 71) and a, (T) = ¢ =a, (T*) by Theorem 10. 


Theorem 12—If T is an iso-abelian operator on a smooth reflexive Banach space 
X, then ¢ (T) 1 V(T) C op (T). 


Proor : It is clear that o (7*) lies on the unit circle. Since X¥ is smooth and 
reflexive, X* is strictly convex, and therefore by (Bonsall and Duncan’, Th. 8, p. 93) 
(T*) VV (T*) C 5, (T*). Since V(T) C V (T*) and o (T) = 3 (T*), ¢ (T) QV (T) 
C o(T*) (1) V(T*). The proof is completed by Corollary 11. 


Corollary 13—If T is an iso-abelian operator on a smooth reflexive Banach space 
X, then the extreme points of Co o (T) which lie in V (T) are the eigen values of T. 


Theorem 14—If T is a normal and iso-abelian operator on a smooth reflexive 
space X, then V (T) is closed if and only if the spectrum of T consists entirely of the 
point spectrum. 

Proor : Since T is an invertible isometry, « (7) lies on the unit circle. As T is 
normal, V (T) = Coo (T) and therefore it follows that E = «(T). Suppose V (T) is 
closed. Then by Theorem 1, E 1 5. (T) = ¢. But this holds only if « (7) = o, (T) 
because o, (7) = ¢, by Theorem 10. 


The converse is an immediate consequence of the fact that T is convexoid and 
Co o,(T) C V (T) (Bonsall and Duncan’, Th. 3, p. 21). 


Theorem 15—If T is an iso-abelian (or normal) operator on a Banach space with 
descent @(T — +/) finite for each » in « (T) then Co co (T) C V (T). 


Proor : By Corollary 3.6 of Mattila’, the ascent « (T — J) < 1 for each pin C 
and hence @(T — pl) = «(T — pl) <1 for each »pino(T). Now pz is an isolated 
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point of o (7) and is a pole of the resolvent of T (Lay’, Th. 2.1). Therefore p» is an 
eigen value of T (Taylor’’, Th. 5.8 — A). Hence s(T) = op (T). Since Co op (T) 
< V(T) (Bonsall and Duncan’, Th. 3), the proof is completed. 


Corollary 16—Under the hypothesis of the Theorem, o (T) C W (7). 


Corollary 17—If T is a normal operator on a Banach space with @ (T--p/) finite 
for each » ins (T), then V (T) is closed and convex. 
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Yadav’ proved ¢, (t) € BV (0, =) implies 3 An en E | N, Pn | under certain 
conditions on {pp} provided {en} € BV (0,0). In the present note, we 
have extended the above result to matrix summability. 


§l. Let 2 u, be a given infinite series with sequence of partial sums isa Lee 
||7|| = (4,,.) be an infinite triangular matrix with real constants. Then sequence-to- 
sequence transformation 


n 


= x Anyk Sky 1 = 0, I; 2; tee al ian 
k=0 


defines the 7-transform of the sequence {s,}. The series 2 u, is said to be T-summable 


to s, if lim ¢, = s, and is said to be absolutely T-summable or simply summable | T | 
noo 


bead . 
if the infinite series 2 | ¢, — t,_, | is convergent. 
n=) 


If matrix element a,,, = 0, for every k > n, then the matrix is called triangular 
matrix. The matrix 7-reduces to Nérlund matrix generated by the sequence of coeffi- 


cients {p,}, if 


p> ifk <n; 


On if keane 


Anysk = 
where 
P= >» Pe F 9. 
k=0 


§2. Let f(t) be a function integrable (L) over (—7, 7) and periodic with period 
27. We assume, without any loss of generality that the constant terms in the Fourier 


series of f(t) is zero, so that 


{ f(t) dt =0 
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and 


F (Gia (adGount.-+ by sin nt) = Ae: (2.1) 
n=1 a) 


We write 
dH=—4t(fe+e+f — Ol] 
Dy (t) = ¢ (t); 


o.() = 5 > | 1 ye Gab dae ees 


0 


da (t) = I' (a + 1) t-* Da (t), a > O; 





yt 
D, (1) = 4 + cost + cos 2+... + cos mt = SRC 4 Dt 


> —~k+1 
Ansk a Gnor: Vary = ARE 
nok 


=F 


Ax QAnysk = Ansk — Ansk+1 
« = ([1/t], 
where [4] is an integral part of A. 


Furthermore, A will denote an absolute constant and we will write A + A = A 
and A —- A= A. 


§3. Bhatt’ proved the following : 


Theorem A—Let {p,} be a non-negative and non-decreasing sequence of real 
numbers such that 


{Pn41 — Pn} is ultimately monotonic, (34) 
] 
— a c BV ...(3.2) 
Py, - I 
bs > P, <C Ki(3.3) 
n=k 


for kK = 1, 2, ..., where C is a fixed positive constant. If ¢, (t) € BV (0, 7), then the 
series (2.1), at f = x, is summable | N, Doin 


The special case of this theorem in which (N, P,) is (C, «) with « > 1, is included 
in a more general result proved much earlier by Bosanquet?. 


ABSOLUTE MATRIX SUMMABILITY 185 


Recently Yadav’ generalized the above theorem as follows : 


Theorem B—Let {p,} be non-negative and non-decreasing and satisfy (3.1), (3.2), 
(3.3) and {e,} € BV (0, ©9), then 


$i (t) © BV (0, 2) > E e, An (t) © | Ns Pa | - 


The problem of absolute matrix summability of series (2.1) where ¢ (t) € BV 
(0, 7) has been studied by Kishore and Hotta‘, Varshney®, Kuttner and Sahney® and 
Jurkat et al.® 


In this note, we prove : 


Theorem 1—Let ||T|| = (@,,«) be an infinite triangular matrix with a,,, > 0; 
Ans0 — 1 Yn > 0 such that 


{Ansk}pno ...(3.4) 
is non-negative and non-increasing sequence with respect to k 
{Anyk 73 ra pee eR KTR) 


is either monotonic non-increasing or non-decreasing sequence with respect to k; 


r—l 
ES sup. | Ac Ven | = O() 3.6) 
k=l k<nSr 
mcd ‘ 
Ansn-k & C (k + 1), for every k aac) 
n=k 
and 
An—15k 2 (1 “F nso) Ansk+1 Bvto.5) 


for every0 <k<n,n> 1, where C is a fixed constant. If ¢, (t) € BV (0, =), then 


the series ¥ e, A, (t), at t = x, is summable aed where {e,} € BV (0, ©9). 
1 


We note that in the case of Norlund summability (3.8) is automatically satisfied 
while condition (3.4) implies that the sequence {p,,} is non-negative and non-decreasing. 
Also condition (3.5), (3.6) and (3.7) reduce to condition (3.1), (3.2) and (3.3) respectively. 


We shall require the following Lemmas : 
§4. Lemma® 1— Let ||7|! = (a,,«) be an infinite triangular matrix with a,,. > 0, 
Ano = | and let {Anz} also satisfy (3.6), (3.8) and 


Bava (4.1) 
rat fe -— C. 


nak 
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If 2% u, is bounded then a necessary and sufficient condition for i to be summable | T | , 
is that 


» weal > K dyyk Ue| < 09. 


Lemma 2—If {b,,«} is non-negative and non-decreasing sequence with respect to 
k, then forO €a<d< 0,0 <t < mand for every n anda 


d 
| 2 47 Da pak eft | x APsipes 
where 
Bysk = % byyr and {ex} € BV (0, 0). 


PROOF: We have e, = «,% — B, such that {xx} and {8%} are non-negative and non- 
increasing for {ex} € BV (0, oo), 


Now 


d 
| > ex Ls et(n—k)t | 
a 


d d 
= | P aK | eae et(n—k)t sence Ss Bx Day et(n—k)t | 
a 
a 
S%_ max | = b,,,_~ etln-*ye | 
aSd\<d a 
=i 
v Ba max | 2 Nea hile lal h | 
a<d\<d a 


S Ai Bisn=s <r AsBasn_t =A Bases 


by using a Lemma‘. 


Lemma 3—If {4,,x};_. iS a non-negative Sequence such that 4,,, = 1 and satisfy 
(3.6), there exists a constant such that, for all 
Var SMO SK <n. (4.2) 


PRrooF : If a,,, > 0, Ano = | and satisfy (3.6) then for 0 < k <n, we have 
Vaaae < M. 


ABSOLUTE MATRIX SUMMABILITY 187 
Since 


Ve, V, = — 2Gnyn—1__ = 


n= lyn 
Qnsn-1 i non 


the result hold fork = 0,1. If2 < k <n, then 


k-1 
| Veseok a eG n-1 | < > | A, V weiter | 


Lac | 


ee Apveenl nO (i) 
fa) 
in view of (3.6). 


We note that, taking the special case kK = n, we have 


M 
Serer TL (4.3) 


§5. Proof of Theorem 1—We know’ that if ¢, (t) € BV (0, 7), then the series 
(2.1) is convergent. We observe that (3.7) implies (4.1). In fact 





nso < 


co 


co 
Axtk 1 >) 
wiih 21.22 ———— ok sa Ge 
n+ 1 = Anon-k G 


n=k 


n= 


Thus by Lemma ||, it is sufficient to prove that 


>: AA oe (5.1) 
n ] 


| 


4. 

where 

6, = > k Ck Ansk Ax (x). 

k=) 
Now 
2 
Ax (x) = s | ¢ (t) cos kt dt 
0 
car ae =| sin kt d di (t) + =| t cos kt dd (t) 
7 k ae 
0 


by integration by parts. Hence to prove (5.1), it is enough to show that the series 


us 


>! GED | | dd, (t) > ex a,x sin kt | ...(5.2) 


am 
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and 





5 = + i | [16 o> Kk ek Gy,x cos kt | (5.3) 


n e k=1 


are convergent. But by hypothesis 
if 
J | dg, (t) | <A 


it suffices for our purpose to show that uniformly forO0 <t<7 


co n 
>) 1 > ; 
i, = acto | €k Gy,k SIN kt | Ch 


om | k=l 
and 
co 1 n 
l= Sa | D> bee ane cos kt | < C. 
n=) k=1 
Now 
ae oe 
nw a> 
1 = 
n< k=0 
l 2 ; 
+ me n+1 | bab €k Qnyk SIN kt | 
ios = 
<Saas 
is ero ad Onsk 
net k=o * 
+ Ax ae : 
n+1 7° max | oars 
n>- O<r<n sin Kt | 
k=O 
l 
= f Az O —______. 
> a o> (n +1) 
v n>t 


O (1) 


in view of (3.6), (4.3) and Abel’s Lemma. 
This completes the proof of (oa). 
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We now proceed to prove (5.3). Using Abel’s transformation, we have 


~ 
hat > oy | SS K ey Ay,x cos kt | 





Sa | k=) 
oe) n k 
1 
sles ‘> | > (rer ane) C> r cos rt) | 
n=1 k=) r=1 
= » + 2 
nwt n>¢t 


= S, (t) + S, (t), say. 


Using (3.4), we observe that 


S, (t) t <)>: Tri > | Ax (€k Gnsk) | 


n=) 


< t PE n 2 | (ex Ax Qnsk 9 Anyk +1 A ex) 
=1 


n=l 


tr n 
it a |e Ag Byte 
k™) 


neal 


+tin >» | Anka A ex | 
1 k=1 
tr n r n 
gKAt Sn [Aran | tt ona, = |Ae| 
n=1 k=1 1 k=l 


r n i 
= At > nS (nk Pac Anyk-+1) a At % A Ans 
k=0 


n=1 


r 
< At 2 Tt Gag 


0 (1): ...(5.4) 


Now since 


5 rcosrt = O(t-?) +nD,(t)0<t<7 
1 


r= 


we have 


n k 
S,(t) =t 3 aren} | S, [Ak (€k Gnsk)] > rcor rt) | 


n>t kw Se 
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; n 
oo O[rt > ain > | Ax (€x Gn k) | 


n>t k 


+ t >a | > Ax (ex Gnyx) k Dy (t) | 
n> 1 
= O [rt > any», Ak (ex Onsk) 
n>t k=l 
+ L> ea | > €k (Ak Gnjk) K De (t) | 
n>Tt 1 
1 n 
+ > ary aT} | > Qnyk+1 (Aex) k Dy, (t) | ] 
n>T a 
= K, + K, + Ks, say. AO} 
Using (3.6), we get 
1 n 
Ky = Ole STE SY | belercnn) 13 
n>t tae | 
] 
=O Dae + 1) Ans 0) 
nu>t 
a Po 
Ol > aaa 
n>Tt 
am (1) (5.6) 


asin S; (t). In order to deal with Ks, we consider two cases separately. 


ae ed he Se, n— 
Case I—Let {Ag any}, = {Bayx} *_ be monotonic non-decreasing with respect 
tok. Then we have by Lemma 2, and Abel’s Lemma 


Ee (4, a,,r) r D, (t) 


= 12 er Bar 7 Dy (t) | 


Thus 
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n—1 
i + 1/2)t 
< | S Cr Daye See ee Let] Ren Ga_ D, (8) 
r=] 


k 
o(+) max | whe sin (r + 1/2) ¢ | 
I1<k<gn—] | — 


FROUS aaa lier'|) 


O (n/t) Pe 3 Gy 55] 
O (n/t) (Gay_—2 ze Bisa) 


n Qnyn—F a’ 
K, = O [e& (eta) Ttaw)] 


O[ x Cpe ae Gn sn) 
n>t 


A, 


= o[ x Se ie Ansel 
n> n>7T 
= O(1) pen) 


in view of (3.7) and (4.2). 


Case II —Let {Ak (4ysk)},29 == {B,,%}22, be monotonic non-increasing with re- 


spect to k. Now since =r D, (t) = O (n/t?), so we see that inner sum in K, is 
[mee | 


Hence 


But 


= O [n 72 (4,0 cin noid]. 


K; =O [t py (4,50 =P Qns1)]. 
n>t 











1 NAnsi NAnyi ke ve 
G,59 — 2an = eral [2 + 1) 4,50 — ae a5 n+l was 1] 
Anyi 
~~ nti 
l NAny RGyy ey Cask oe 
= n ae 1 [(n + 1) Qnyo — Aak ] (n ze DA Anso n+ 1 


; 1 
= EH Mn + Dame — BI + 0 (a »): 
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Thus using (3.6), we have 
1 a ] 
K, = O[t > CRS) [ (2 + Iau — {311+ Olt S n(n-+1) 1! 














A>t n>T 
= 0O(1) ...(5.8) 
so that K, = O (1). 
Also 
ke > st», M(B: einen eDen 
n>t : k=) 
<>; ty > 1 Aves ent 
n>T 
co 1 n 
Sr Flidee (ase 
n=) kel 
co co 
-01S rina Stns 
tom] n=k+1 
But 
2k 2k C 
Ansk+1 
Ss mone Tes YE > OInket S EF 
neak+] n@kt+1 
in view of (3.4) and A,,o = 1. We have, by Lemma 3 
por ie SNE oY 
i. €. 
M Ayyx M 
and so 
acd fore) 
Ansk+1 1 1 
SS < M > air = O= 
ite seh dea pe (n — k) (n+ 1) (k + 1) 
i 6: 
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Thus, we have 
oo 
K, = O A. | ex |) = OX1). (5.9) 


Combining (5.4) to (5.9), we get 
I, = O(1). 


This completes the proof of Theorem 1. 
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SCATTERING OF COMPRSSIONAL WAVES BY A CIRCULAR 
CYLINDER 


P. Kesari AND B. K. RAJHANS 
Department of Applied Sciences, Indian School of Mines, Dhanbad 826004 


(Received 23 March 1985; after revision 9 September 1987) 


This paper deals with the scattering of compressional waves by a circular 
cylinder. The cylinder is embedded in an unbounded isotropic homogeneous 
elastic medium and it is filled with some acoustic fluid.T he line source, gene- 
rating the incident pulse is situated outside the cylinder parallel to its axis. 
The problem is investigated by the method of dual integral transformations. 
The resulting integrals are evaluated asymptotically to obtain the short time 
estimate of the motion near the wave front in the illuminated region of the 
elastic medium. We interpret the approximate solution in terms of geometri- 
cal optics. 


1. INTRODUCTION 


Lapwood‘ considered the problem on the disturbance due to a line source in a 
semi-infinite elastic medium and obtained the exact formal solution of the problem. 
Jeffreys and Lapwood!® discussed the reflection of a pulse within a sphere and obtained 
the solution using operational methods. The problem on Scattering of two dimen- 
sional elastic waves with a cylindrical obstacle in unbounded medium has been consi- 
dered by various autors in recent years. Gilbert and Knopoff? used Friedlander’s® 
method to investigate the problem of scattering of impulsive elastic waves by a rigid 
circular cylinder. Gilbert® discussed the problem of scattering of impulsive elastic 
waves by a smooth convex cylinder using the Same method, Mishra!516 applied Fried- 
lander’s method to investigate the problem of scattering and diffraction of two dimen- 
sional sound pulses by a acoustically semi-transparent circular cylinder. Hwang 
et al."* applied a similar method and discussed the case of three dimensional elastic 
waves scattring by a rigid cylinder in an elastic medium. 


ing elastic medium at a finite distance Parallel to the axis of cylinder, We assume 
that the velocities of P and SV waves outside the cylinder are « and 8 respectively and 
that of P-waves inside the cylinder is %. To be specific, we suppose a > a > Bp. 
This assumption of the velocity distribution corresponds to the actual velocity distri- 
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bution of elastic waves inside the earth and to the location of the source in the mantle 
and the outer core as the obstacle!2. We also suppose the density of the medium 
outside the cylinder is P and that inside the cylinder is P’ where P > Pf’. The present 
discussion therefore may have some relevance in seismological problems. 


2. FORMULATION OF THE PROBLEM AND FoRMAL SOLUTION 
Let the axis of the cylinder be taken as the z-axis and let a co-ordinate (r, @) be 
located in the (x, y) plane with @ = 0, r = ro, (ro >a) corresponding to the location of 
the line source which is parallel to the axis of the cylinder. The equation of the 
cylinder is r = a. 


The governing equations for the present case are 








1 @o - 2” 
oe ore —_ Vv? — 7a 8 (r a ro) $ (t) ) (6), (r > a) (231) 
1 au . .: 
muna Te v7 = 0 (rp 2a) eAZ2) 
ine = 
oar ~ Vm =0, (r <a) -+(2.3) 


where, y? is the Laplacian operator. 


Besides the potentials also satisfy the conditions of quiescence at t= 0. The 
boundary conditions for the problem are that the tangential and normal components 
of stress vanish outside the obstacle and the radial displacement is continuous at the 
boundary of the obstacle'’. 


Applying Laplace and Fourier transformations’ to (2.1), (2.2) and (2.3), and 
following Mishra'*'*, we find that Laplace transforms of the solutions are given by 


co 


9 (r, 0,5) = | Tp ( 7) Ky (22) exp Ciné) a 


ny i Ky ( =) Ku (= lr exp (i16) du, (ry > r > a) 
; (2.4) 
U(r, 0, 8) = ( Ku ( A ) Ki (“2 ) * exp (ip) du, (r >a) ...(2.5) 


and 


co 
Go (r, 8, 5) = | LR + NP Tu (= ) Kp (22) exp (i6) dp, (r < a). 
-00 ...(2.6) 
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Here L, M, N, P, O and R are the values of the constants which are determined 
with the help of boundary conditions. 


We see that (2.4), (2.5) and (2.6) give the integral representation of Laplace 
transform of the formal solutions. The time solution can be obtained on performing 
Laplace inversion. 


3. INCIDENT, REFLECTED AND REFRACTED PULSES 


We first give a brief description of the geometry of the problem. Initially the 
incident P-pulse striking the outer surface of the cylinder gives rise to reflected P, 
reflected S and refracted P-pulses, acording to the laws of ordinary geometrical optics 
(Fig. 1) when the rays strike the outer surface of the cylinder at critical angle, the 
reflected P-rays become tangential to the surface and as such they move along the 
surface. These surface waves at each point of their path give rise to SV-wave at 





Fio. 1. (1) Incident rays, (2) reflected rays, (3) refracted rays and (4) diffracted rays. 
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critical angle in the outer medium and P-wave along the tangent to the surface". The 
former are denoted by P (P) S and latter by P (P) P (Gilbert and Knopoff’). In the 
wie of grazing incidence, the disturbances move along the surface and at each point of 
their path they shed diffracted P-waves in the outer medium trangential to the surface. 
These are denoted by P (P) P (Bullen'*). The reflected pulses are denoted by P. P and 
P-S and refracted pulse is denoted by PKP (Bullen?*). 


Now, we use the saddle point method? to obtain the solution in the illuminated 
region of the elastic medium. Therefore using the various approximations for modi- 


fied Bessel functions as given by Mishra'®* to the integral (2.4), (2.5) and (2.6) we 
find 


ae or 1/2 a \1l2 
@ (r, 8,5) & ( ) exp (—st,) + (= ) A, B, exp (—Stz) 


(ro >r Da) ...(3.1) 
y (r, 0, 5) (= Ne A, B, exp(—sts), (r > a) e432) 
and 
ot 1/2 
Go (r, 0, 5) (=) Az Bs exp (—St,), (r <a) (3.0) 


where, ¢,, t,, t3 and ft, are respectively the arrival times of the incident, reflected and 
refracted pulses. 


4pe 


52 Sin® y cos 1 (1 — nj sin? nj? (1 — n? sin? »)'? 





+ to P’ cos  — Pa (1 — 2n? sin? »)* (1 — ni sin? y)'/? 


bE = 4p B*a-2 sin® y cos n (1 — nj sin’n)'/? (1 — n?® sin? y)1/2 


+ P’ a cos y + Pa (1 — 2 n® sin® y)? (1 — n? sin? n)1/? 


ae aa COs 1 ‘3 
1” Lo Rs cos & + r Ro cos f 


A,, B, and A;, Bs have similar expressions as above. 


R,, R2, Rs denote distances between the source and receiver and the angles are 
defined in Figs. 2, 3, 4 and 5. 


Now we can obtain the short-time approximations for the solution by perform- 
ing Laplace inversion“. We find that 





H Le t; H(t sa ty) ~ 
9 (r, 6, t) ad 2t, is a + Ai B, (t — hp? (ro 2T2 a) rae) 
bile, 0; 9).00 Ay By 2D, Or & a) 3.5) 


(t Ln tz)1/? , 
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Fic. 2, Geometrical interpretation of the saddle Fic. 3. Geometrical interpretation of the 
point for the incident P pulse. saddle point for the reflected P pulse. 





Fic. 4. Geometrical interpretation of the saddle Fic. 5. Geometrical interpretation of the 
point for the reflected S pulse. saddle point for the refracted P pulse. 
H(t — t,) 


90 (r, 6, t) © As Bs 7? (r <a) (3.6) 


Figures 6-8 present the numerical evaluation of these pulses respectively. 


COMPARISON 


Hwang et al.’* discussed the problem of three dimensional elastic wave scattering 
due to a rigid cylinder in case of a compressional point source on comparing our 
results with those obtained by them, we find that in addition to the results obtained 
by us, they obtain an additional event in the illuminated region which they term as 
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Reflection Co-efficient (P-waves) 





Fic. 6. Reflection co-efficients for P waves. 











Reflection Co-efficrent (S-waves) 


O° 10° 20° 30° 40° 50° 60° 70 80° 90 


Fic. 7. Reflection co-efficients for S waves. 


diffracted PPaS wave. It may be pointed out that this possibility does not hold in 
the illuminated region. The physical ground for this is that diffracted waves exist in 
the shadow region only. Diffracted field isin fact the contribution of poles. It is 
obtained on evaluating the integrals (2.4) and (2.5) by Watsons residue method. 
Besides the results obtained by the Saddlepoint method are interpretable only in the 
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2-2 
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Refraction Co-efficient (P-waves ) 
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-0.3 


Fic. 8. Transmission co-efficient for P waves. 


illuminated region’. This conclusion seems quite reasonable. It agrees with the results 
obtained by Mishra!®!6 and Rajhans and Agarwal?. 
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HALL EFFECTS ON THERMOSOLUTAL INSTABILITY OF A PLASMA 


R. C. SHARMA AND NEELA RANI 
Department of Mathematics, Himachal Pradesh University, Shimla 171005 


(Received 6 November 1986) 


The thermosolutal instability of a plasma in the presence of a uniform horizo- 
ntal magnetic field is studied to include the effects of Hall current. When the 
instability sets in as stationary convection, the Hall currents and the stable 
solute gradient are found to have destabilizing and stabilizing effects 
respectively. The case of overstability is also considered wherein the 
necessary conditions for the existence of overstability are derived. 


1. INTRODUCTION 


The theory of thermal convection of fluid layer heated from below under vary- 
ing assumptions of hydromagnetics has been treated in detail by Chandrasekhar’. Gupta? 
has studied the problem of thermal instability in the presence of Hall current. The 
thermal instability of a compressible Hall plasma has been considered by Sharma‘. 
Veronis® has studied the problem of thermohaline convection in a layer of fluid heated 
from below and subjected toa stable salinity gradient, whereas the problem of 
thermohaline convection in a horizontal layer of viscous fluid heated from below and 
salted from above has been considered by Nield*. For thermohaline convection, 
buoyancy forces can arise not only from density differences due to variations in temper- 
ature, but also from those due to variations solute concentration. 


The physics is quite similar in the stellar case in that helium acts like salt in rais- 
ing the density and in diffusing more slowly than heat. The conditions under which 
convective motions are important in stellar atmospheres are usually far removed from 
consideration of single component fluid and rigid boundaries and therefore, it is desi- 
rable to consider a fluid acted on by a solute gradient and free boundaries. The problem 
of the onset of thermal instability in the presence of a solute gradient is of great im- 
portance because of its application to atmospheric physics and astrophysics, especially 
in the case of the ionosphere and the outer layers of the solar atmosphere. The Hall 
effects are likely to be important in these regions. A reconsideration of the thermoso- 
lutal instability including Hall effects is certainly called for and is the object of the 
present paper. 


2. FORMULATION OF THE PROBLEM AND DISPERSION RELATION 


| Here we study the effect of Hall current on the thermosolutal instability of a plasma 
in the form of an infinite horizontal layer of thickness d, acted on by a horizontal 
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magnetic field H (H, 0,0) anda gravity force g (0,0, — g). This layer is heated 
from below and subjected to a stable solute gradient so that the temperatures and 
concentrations at the bottom surface z = 0 are 7) and C, and at the upper surface 
z = dandT, and C,. Consider the cartesian coordinates (x,y,z) with the origin 
on the lower boundary z = 0 and the z-axis perpendicular to it along the vertical. 
The layer is heated and soluted from below such that a uniform temperature 
gradient 8 (= | dT/dz | ) and unitorm solute gradient 8’ (= | dC/dz | ) are maintained. 


a 

Let g (u, v, ¥’), p, P, T, C, «, «', g, n, N and e denote respectively the velocity, 
pressure, density, temperature, solute concentration, thermal coefficient of expansion, 
an analogous solvent coefficient of expansion, gravitational acceleration resistivity, 
number density and the charge of an electron. Then the equations expressing the 
conservation of momentum, mass, temperature, solute mass concentration and 
equation of state are 


> 





ee +9) 9 = - x wtovg+e(l + | 

+ an (Vv x H) x H 2th) 
9 = 0 Q) 
Sf + GN) T= «VT m2) 
eG WC He VEC (4) 
p =p, (1 —a(T — To) + «’ (C — Cy) ...(5) 


where the suffix zero refers to values at the reference level z = 0 andin writing eqn. 
(1), use has been made of the Boussinesq approximation. The magnetic permeability 
jie, the kinematic viscosity », the thermal diffusivity « and the solute diffusivity «’ are all 


assumed to be constant. 


Maxwell’s equations give 


2H Gig) g+nvH— ye VX UW x HD x HI ...(6) 
v.H = 0 alg? 
where _ = — + 0.9 stands for the convective derivative. 
The steady state solution is 
9 = 0,1 = To—Bz,C Cy — Bz, P = p,(1 + 282 — «’B’z). —---(8) 
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Consider a small perturbation on the steady state solution and let q (u, v, w), 5p, 


=? . 
5p, 6, y andh (hx, hy, hz) denote respectively the perturbations in velocity, pressure, 
density, temperature, concentration and magnetic field. The change in density 5P, caused 
by the perturbations @ and y in temperature and concentration, is given by 


5p = — Po (ad — a’ y). ...(9) 


Then the linearized hydromagnetic perturbation equations become 








1 = — yap + wy'g — 8 (09 — «'y) + inp (VY XA) x HL ...(10) 
v. q=0 (11) 
oh sation 0 = I] Looe Ey 

Fre SLI YC Ble A dnNe VX (VY x A) x H] (12) 
v.h =0 fell d) 
2 = Bw + « 96 ...(14) 
Y= Bw te yty, | (15) 


Analyzing the disturbances in terms of normal modes, we assume that the per- 
turbation quantities are of the form 


[w, 0, y, hs, &, &] = [W (z), © (z), F(z), K (2), Z (2), X (z)] 
exp (ik,x + ikyy + nt) st LG) 


where n is the growth rate; kx, ky are wave numbers along the x and y directions 
and k = («: + k* es is the resultant wave number € = wad and & = ats 
ox oy Ox 


Ohx : 
a naa stand respectively for the z-components of vorticity and current density. 


Letting a=kd,o = nd*/v, Pi = |x, Ps vin, q=y', D= - and putting 
: Zz 
the coordinates x, y, z in the new unit of length d, eqns. (10) — (15), using expres- 


sion (16), in nondimensional form give 


22 
(D* — a*) (D? — a? — o) W— (224 ) @o —aT) + ik, eftd* 
4nPov 
x(D? — a’) K = 0 (17) 
He Had? 


(DS ate= o) Fete Hd* 
°) te ei ...(18) 
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Hd" Hd? 
D? —a* — K =— ik, j 
( P2 S) ik ga ye x ...(19) 
ete f be : Had? . 
(D2 — a? — p,c) X = — ikz 7 Z — ik, 4aNe, (PO) K 
...(20) 
2 
(Dt — q*\—'p.0)'8 == — ik Ww on 049) 
442 
(D? — a? — go) T =— Be W. area) 





Eliminating Z, K, X, © and I between eqns. (17)—(22), we obtain 
(D? — a?) (D*? — a? — oc)? (D* — a? — po)? (D? — @ — ps) (D*? — @ 
— qs) W + Q a*cos* 9 (D* — a*) (D? — a? —p,c) (D2 —a? — qs) 
{2 (D? — a? — oc). (D? — a? — pc) + Oa? cos? 5} W—M a 
cos? 9 (D® —a?)? (D? — a® — a)? (D® — a® — ps) (D? — a? — qo) 
W + {Ra® (D? — a® — qs) — Sa® (D? — a® — p,o)} 
[(D? — a® — pc)? (D*? — a? — c) + O a? cos? 6 (D? — a? 


— ps) —M a’ cos? 9 (D* — a*)(D? — a? — s)|W=0 
...(23) 


where M = (H/4rNen)* is nondimensional number accounting for the Hall currents, 


Q = ue H* d*/47Povy is the Chandrasekhar number, R = gaBd*/v«, is the thermal Ray- 
'Q' 4 


ga’B'd 


leigh number, S = ats is the analogous solute Rayleigh number and cos 9 = k;/k. 


We consider the case of a fluid layer with two free surfaces when the adjoining 
medium is assumed to be electrically nonconducting. The case of two free boundaries 
is the most appropriate for stellar atmospheres’. The boundary conditions appropriate 
to the problem, using (16), are 


W=DW=X=DZ=-0=Lf=0 ...(24) 


> 
and the components of 4 are continuous. Because the components of the magnetic 
field are continuous and the tangential components are zero outside the fluid then 


DK = 0 A) 


on the boundaries. 


Using the boundary conditions (24) and (25), it can be shown with the help of 
eqns. (17) — (22) that all the even order derivatives of W vanish at the boundaries. 
Hence the proper solution of eqn. (23) characterizing the lowest mode is 
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where W, is a constant. 


Substituting (26) in eqn. (23) and letting x = @’/n®, io, =o/n°, R, = Rin", 
S; = S/n* and Q, = Q/x?, we obtain the dispersion relation 


(1 + x) (1 + x + ipo) (1 + x + ios) (1 + x + ip.) 


+ Q,x cos? 6}? + Mx cos? (1 + x) (1 + x + is;)*] 
+d + x + ip, 1)? (1 + x + io,) + Q, x cos*a(1 + x 


+ ip, 6) + Mx cos? 6(1 + x) (Hx + i,)] 


R 


S; (1 + x + ip, 6) 


re (1 + x + igs,) 


(27) 


3. THe STATIONARY CONVECTION 


For the stationary convection, ¢ = 0 and eqn. (27) reduces to 


1+ x\ {(l + x)? + Q: x cos? 6} + Mx cos?6 (1 + x)? 
Rea oy | =—— fT oe x lc coat aes cote (ieee 
(1 + x)* + Q, x cos? 6 + Mx cos? 6 (1 + x) 

.-.(28) 

Equation (28) expresses the modified Rayleigh number Ri asa function of the 

dimensionless wave number x and the parameters S,,Q, and M. To investigate the 

effect of stable solute gradient and Hall currents, we examine the nature of dR,/dS; and 
dR,/dM analytically. Equation (28) yields 


aRee 
ie ...(29) 


which implies that the stable solute gradient has a stabilizing effect on the thermoso- 
lutal convection. Equation (28) also yields 





dR, = — (LEX) BC a cost HC 4 a+ 01 x cos" dl 
dM o (1 + x)? + Q, x cos? 6+ Mx (1 + x) cos? 6}? 

: ...(30) 
which is negative. The Hall currents, therefore, have a destabilizing effect on the 
thermosolutal instability of a plasma. 


4. THE OVERSTABLE CASE 


Here we discuss the possibility as to whether instability may occur as _ overstabi- 
lity. Since for overstability we wish to determine the critical Rayleigh number for the 
onset of instability via a state of pure oscillations, it will suffice to find conditions for 
which (27) will admit of solutions with o, real. Equating real and imaginary parts of 
eqn. (27) and eliminating R, between them, we obtain 


A,c* + Asc® + Aye? + Ac + A, = 0 eal Ga 
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where we have putc = oa; , b = 1 + x and 


A, 
Ag 


piq’ (1 + pi) 6 ...(32) 
(1 + pi) bo + 2M (1 + pi) (b — 1) BF cos? 6 
+ b7 {Q, cos? 6 (b — 1) (3p1 — P2 + 2) + M? cos* @ (1 + Pi) 
x (b — 1)?} + b8 {S, (p, — 9) (6 — 1) + MQ, cos* 6 (3p; 
+ p, + 2)(b — 1)?} + 5° {Qi cos* 6 (3p1 — 2P2 + 1) (b—1)? 
+ 2M S, cos? 6(p, — 9) (b — 1)?} + 5 (20, Si cos” 6 
(p1 — p2) (b — 1)? + M?S, cos‘ 6 (p,; — 9) (b — 1)’ 
+ QO? Mcos*6(pi — 1)(b — 19%} 
+ b® {Q3 cos® 6 (pi — p2)(b — 1)? + 2M Si Q, cos* 6 
(p: — 9) (b — 1°} + QF Si cost a (p, — 9) 5 (b — 1%. (33) 


Since «, is real overstability, the four values of ¢ are positive. The product of 
the roots of (31) is A»/Aa, which is positive if Ay > 0 (since from (32), As > 0), 
Equation (33) shows that Ao is always positive if 


I 


Pi > 1, Pp, > qand pi > Pz apt o4) 
which means 
k<yv,« <« and« <q s(32) 


Thus « < v,« < «’ and« < are the necessary conditions for the existence of 
overstability 
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The non-linear partial differential equations of motion for three-dimensional 
flow of an incompressible viscous fluid flowing over a semi-infinite plate under 
the influence of a magnetic field anda pressure gradient, and with or without 
suction or injection through the plate wall are developed with necessary 
boundary conditions. The group-theoretic methods are employed to obtain 
proper similarity transformations. Similarity solution of such flow system 
yields coupled non-linear ordinary differential equations. 


1. INTRODUTION 


Due to increasing number of technical applications using magnetohydrodynamic 
effect, it is desirable to extend many of the available viscous hydromagnetic solutions 
to include the effects of magnetic fields for those cases when the viscous fluid is electri- 
cally conducting. Flow past a flat plate has been studied in Rossow! and Carrier and 
Greenspan’. Heat transfer for these cases has been discussed in Rossow? and Afzal?, 
Rossow' has considered transverse magnetic field whereas Carrier and Greenspan? have 
studied the effect of a longitudinal magnetic field on the velocity and the temperature 
distributions. There have been other computations of viscous magnetohydrodynamic 
flow reported in the literature*’®, Djukict has considered only Hiemenz magnetic flow 
of non-Newtonian power-law fluid. He has solved the governing non-linear ordinary 
differential equations using Galerkin’s technique, whereas Srivastava and Usha® has 
discussed magnetofluiddynamics of simple two-dimensional flow with pressure gradient 
and fluid injection. Recently such case is extended by Timol® for a non-Newtonian 
power-law fluid. 


2. Basic EQuaTIONs 


The basic equations of magnetofluiddynamics and coventional fluiddynamics are 
different by only additional force term due to electromagnetic field in momentum 
equation and a term due to Joule heating in the energy equation. In such situation 
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the Maxwell’s equations have to be satisfied in the entire flow field as well as in the 
body and interface. 


In order to derive the basic equations, the following assumptions are made : 


(1) The fluid under consideration is incompressible finitely conducting with 
constant physical properties. 


(2) Hall effect, electrical and polarization effects are neglected. 
(3) The induced magnetic field is neglected. 


(4) The flow is steady and laminar and the imposed magnetic field is perpendi- 
cular to the free stream velocities. 


(5) The magnetic Reynolds number is assumed to be small. 


Under these assumptions we now write continuity and momentum equations 
governing the velocity distribution in the presence of magnetic field as, 


v.V=0 rat i) 


PV. W = — VP + Pyro VV +Ix 8B P84, 


where the third term on the right-hand side of eqn. (2) is the Lorentz force due to the 
magnetic B, and is given by 


J x B=o(V x B) x B. ...(3) 


We now consider the flow past a semi-infinite flat plate placed in the direction of 
the flow. The plate is in ¥ — Z plane and is between 0 <x < coand -w< Z 
< co and free stream is in the X-direction as shown in Fig. 1. Here we shall consider 
that all flow quantities are independent of the Z-coordinate. Such flows are charact- 
erized by the fact that their stream lines form a system of ‘translates’. That is stream 
line pattern can be obtained by translating any particular stream line parallel to the 
leading edge of the surface’. It is hoped that by assuming independence of flow 
quantities in one direction, more quantitative information may be obtained on the 
characteristics of three-dimensional boundary layer flows. 


Thus the problem considered here is essentially a quasi-two-dimensional one. 
The magnetic field vector B is perpendicular to the free stream and is along Y-direction 
because the induced magnetic field is neglected. The basic equations of the steady, 
viscous, laminar, there-dimensional boundary layer flow of a fluid flowing over a semi- 
infinite flat plate under the influence of transverse magnetic field and pressure gradient 
with or without suction or injection through the plate wall can be derived from the 
eqns. (1) and (2) on the basis of usual boundary layer assumptions’, as 


Ou ov 0 - (4) 


ox oy 


210 


Fic. 1. 


M. G. TIMOL AND1A. 


G. TIMOL 


Vv 


u,U 





Flow over a plate in rectangular coordinate system under the influence of transverse 


magnetic field. 


au, ou _ og dp 
ow ow sg: ap 
a GL p dz oy lle 


Ce 9 
Be — [go Be (x)/P] u 


2 
a — [go BY (x)/p] w, 


The boundary conditions for such flow system are taken as, 


Case I 


Case II 


where x and z are the coordinates parallel to the plate, 


y=0>u=—0;0 = Vy x™-/2, yw = Q 


y—aco > u—> Uo x"; w > W (x) = 


2 


ie aah 


- 


y=0>u=0;v=V, exp (™); w =o 


yroos>u >U (x)= Uy exp (mx) 


w —> W(x) = Wo exp (mx) 


aa) 


.. (6) 


PS ah 


. (8) 


...(9) 


...(10) 


ya distance from the plate; 


u, v, W the velocity components of the fluid along x, y and z directions respectively; 


U, W velocity components in the main flow alon 
acceleration of the gravity; 


(= m/p) P 


& x and z direction respectively; g the 
P the pressure; yq the magnetic viscosity of the fluid 


—the fluid density; m a physical constant, o (Greek ‘sigma’) the electrical 
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conductivity of the fluid; By the imposed magnetic induction parallel to y axis; Uy, Vo, 
W,, m > 0 are constants. 


Now introduce stream function v (x, ») such that 


ea) eee ee 

u= i Pe ax LL) 
Then, the equation of continuity (4) gets satisfied identically and equations of motion 
(5) and (6) with boundary conditions (7)-(10) will be, 





os ih i ar =Gi (x)+7H oe — H, (x) i ...(12) 
as Fe TO tn — H, (x) w ...(13) 
where 
Gi (x) = — g/P ap (14) 
G, (x) = — g/p a ...(15) 
H, (x) = BF) (16) 
Then the boundary conditions are, 
Case I 
y= 0+ 3 = 0; = — Vy xim-0P w = 0 Sl Ww 
yroos + U(x) = Uy xm w > W(x) = My x”, ...(18) 
Case II 
y=O0o3 = 0,52 = — eexp("F"*),w = 0 ...(19) 


yro> a — U(x) = U, exp (mx); w > W (x) = Wo exp (mx). 
Y ...(20) 


The goal of reducing partial differential equations (12) and (13) to ordinary 
differential equations with the meaningful boundary conditions (17) -(20) by eke 
technique will cause restrictions to be imposed on the functions G, (x), Ge hag 
H, (x). In order to obtain suitable similarity variables, the group theoretic metho 


is employed. 
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3. Group-THEORETIC ANALYSIS 


Similarity analysis by the group-theoretic technique is based on the concept 
derived from the theory of transformations of groups. Recently this technique is 
found to give adequate treatment of boundary layer equations. The basic concept of 
this method was first introduced by Birkhoff® and later on it was extended exclusively 
by Hansen’, Morgan", Ames’ and others. For the present problem we shall employ 
two distinct classes of one parameter groups of transformation—a linear group of trans- 
formation and a spiral group of transformation and also we shall discuss both the 
cases separately. 


Case I] 
A one parameter group of transformation I; is selected as, 


a, ieee y= A%5 
‘ v A*3 0 w= A%4 0 


where «1, a2, a3, a, and A are real arbitrary constants. 


We now seek relations among «’s such that the basic equations (12), (13) along 
with the boundary conditions (17) and (18) will be invariant under this group of trans- 
formation l’,. This suggests that G, (x), Gs (x) and H, (x) are to be selected as 


G; (x) = Bo x*, ele) 
Hy, (x) = hy x weal 22) 
G, (x) = g, x%. ...(23) 


Now invariant conditions demand that the power of A in each term of the transformed 
equations should be equal, which yields 
2a3 — 2%, — a, = asa; = a3 — 3a, = ag % + as — a, ...(24) 
My by — he — O = ay @ = a, — 2ao = ag ay + Oy, aed ke 
Further, from the boundary conditions (17), (18), we get, 


m ] - ( 
a3/a, = mat. Oy [a1 = et af =m me el) 


Oy 
and from (24), (25) and (26), we get, 
a = 2m — lj ag = m — 1; a = 2m — 1. as(at) 
Now invariants of the group I’; are, 
ay]x*s!*1; b/bx*,/*1 and wiex®,I*, 


where a, b and c are constants. This enable us to derive 


\ the following similarity 
independent and dependent variables : 
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n= ay/x(m-1)/2 (28) 
fi (a) = ¥]b xtmtny12 ...(29) 
fo (n) = wie x™. ...(30) 

Substituting eqns. (28), (29) and (30) in eqns. (12) and (1+) and simplifying, we get, 

= * yHa* hy ‘ 
A= aa apa es Y UF WU, --»(31) 
AREA pee ae aE (32) 
1 lJ? =r sy oT = 5 ap Wy We Ge 2. eee 
where 
2m m+1 
A = at A, = 5 (30) 
Grix}i= 2,32" -0(34) 
G, (x) = g, x°"™") =. (35) 
HS) = he ee ...(36) 


and a prime denotes differentiation with respect to . 


Also from the boundary condition (18), it follows that 








U, = ab, W, = c. ees 

Setting 
ee ee = & ee hy 38 

So ca re U2 » 93 Pagal A» Uy W, and Ho Av U, et ) 
and 

fae ...(39) 

TG ( 
eqns. (31) and (32) will become 

Tardis fit So 1s Jy. = 9 .-.(40) 

t; +hS, —A fi fo +S — profs = 9. (41) 


Now when free stream velocites are function of x, then pressure gradients will be 
given by : 


dp dU , go By (x) w+0(42 
<= TC) pica oe ne rer vamos ie Ll (42) 
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dp _ yw , so By (X) wy 4 ayyay, ...(43) 
— Cl) a te eH 
But for steady flow, free stream velocities outside the boundary layer is independent 
of time i.e. 
ee ...(44) 
ot ee 
under (14)—(16), (17) and (34)—(35), the relationship of (42)—(43) will be, 
So — A, + Ho ...(45) 
S,; = A, + po. ...(46) 


Substituting (45) and (46) in the equations (40) and (41) respectively, we get, 


fi +A +40 -—f2)+ m1 —f’) =0 ...(47) 


f2 tht (l-fi fy) + ve Ul — f) = 0. (48) 


Equations (47), (48) are coupled non-linear ordinary differential equations for the 
function of f(y). With the transformed boundary conditions, 


? —V. 
7=0= Ae “n Br os = 0 ..(49) 


Varo f= 1 f/f = 1. -..(50) 


From equations (37) and (39) values of a, b and c will be 


+ fA, Uy yur YH Uy \2 
ane (“* ) A ge (=) ea (51) 


Case I] 


One parameter spiral group of transformation I’, can be chosen in the form, 
=e seer gan Or 
: v= hob wo Brg 


where 1, B,, Bs, B, and 5 are real constants andeisa real parameter of the group 
transformation I. 


In this case I’, enables us that G, (x), Gz (x) and H, (x) are to be selected in such a 
way, 


G, (x) = g, ef* ...(52) 
G, (x) = gi Pe ...(53) 
HY, (x) = hy e®o* ...(54) 
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following the same procedure as in case-I, the followining absolute invariants are 
obtained : 


E = a, Ylexp (4 mx) ...(55) 
b = b, exp (mx) F; (8) ..(56) 
w = c, exp (mx) F, (&) ARS 
where 
1/2 1/2 
a, = (fs ) eS (72s) Cc, = Wo. ...(58) 


Under (55)—(57), finally eqns. (12)—(13) with the boundary conditions (18), (19) will 
be transformed to following non-linear ordinary differential equation : 


Fv +F, Fe + 20 - F%)4+ (1 - Fi) =0 (59) 
Fi + FF, +2(1-F, F)+ (1 —- %) = 0 ...(60) 


where pty = 2 and prime denotes differentiation with respect to & 
0 


with the boundary conditions 
=2V; 


a . yh ee ee => ...(61 
a tases (a At a Cy 


Fe 00r Som al Py = 1: ...(62) 


Equations (59), (60) with the boundary conditions (61), (62) constitute a pair of non- 
linear ordinary coupled differential equations. 


4. CONCLUSION 
The analysis of laminar, incompressible, three-dimensional magnetofluiddynamic 
boundary layer equations of viscous fluids with stream lines forming a system of 
‘translates’ lead to the following conclusions : 


(i) If we take V, = Oin eqns. (49) and (61), we have magnetofluiddynamic 
boundary layer problem without suction or injection through the plate wall which 
either in Case-I or in Case-II is reduced to a solution of a boundary value problem 
of following third-order non-linear coupled ordinary differential equations : 


fe of foe (bef, ) +e) —f,) 0 ...(63) 


fs +hife" +A, (1 — fq fr) + po (1 — f;) =0 ...(64) 
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where 
A, > 0; po > O. 


The boundary conditions are 


f, 0) = 9; f,' (0) = 0; f, 0) = 0 ...(65) 


fi (ec) = 1, fp (ee) = 1. ...(66) 


(ii) For the non-magnetic case i.e. for Ho = 0, the above set of eqns. (63)—(64) 
will reduce to the set of equations for simple three-dimensional viscous incompressible 
flows obtained by Timol et al.!* (for power index n = | and A, = 1/3). 


(iii) For m = 1, U(x) = U) x, W(x) = W,x, the entire problem will reduce 
to the particular case of three-dimensional magnetic flow near a stagnation point i.e. 
Hiemenz three-dimensional magnetic flows. In such situation the transverse magnetic 
field will be found constant through the plate wall. The present problem is recently 
extended to include non-Newtonian fluids of different models by Timol**, 
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